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Abstract 



In this report we give a detailed account on Hanbury Brown/Twiss (HBT) 
particle interferometric methods for relativistic heavy-ion collisions. These 
exploit identical two-particle correlations to gain access to the space-time ge- 
ometry and dynamics of the final freeze-out stage. The connection between 
the measured correlations in momentum space and the phase-space structure 
of the particle emitter is established, both with and without final state inter- 
actions. Suitable Gaussian parametrizations for the two-particle correlation 
function are derived and the physical interpretation of their parameters is 
explained. After reviewing various model studies, we show how a combined 
analysis of single- and two-particle spectra allows to reconstruct the final 
state of relativistic heavy-ion collisions. 
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Chapter 1 
Introduction 



By now, a large collection of experimental data exists from the first rela- 
tivistic collisions between truly heavy ions [ |135|| , using the 11 GeV/nucleon 
gold beams from the Brookhaven AGS and the 160 GeV/nucleon lead beams 
from the CERN SPS. The first relativistic heavy-ion collider RHIC at BNL 
will soon start taking data at ^/s = 200 A GeV, and in the next decade, 
the already approved LHC program at CERN will explore relativistic heavy- 
ion collisions at even higher energies {^/s = 5.5 A TeV). The aim of this 
large scale experimental effort is to investigate the equilibration processes of 
hadronic matter and to test in this way the hadronic partition function at ex- 
treme energy densities and temperatures. Especially, one expects under suf- 
ficiently extreme conditions the transition to a new state of hadronic matter, 
the quark gluon plasma (QGP) in which the physical degrees of freedom of 
equilibration processes are partonic rather than hadronic [ |151| , |66| , |139| , |109| . 
QCD lattice simulations predict this transition to occur at a temperature of 
approximately 150 MeV [ The experimental confirmation of a possibly 
created QGP is, however, difficult, since only very few particle species, mainly 
leptons, can provide direct information about the initial partonic stage of the 
collision. The much more abundant hadrons are substantially affected by sec- 
ondary interactions and decouple from the collision region only during the 
final 'freeze-out' stage. A successful dynamical model of relativistic heavy-ion 
collisions should finally explain all these different observables, their depen- 
dence on the incident energy, impact parameter, and atomic number of the 
projectile and target nuclei. 

At the present stage, theoretical efforts concentrate on discriminating be- 
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tween different models by comparing them with characteristic observables [ 
|135| , |136| , |76| , . The observed enhancement of strange hadron and low-mass 
dilepton yields and the measured J / \l/-suppression provide strong indications 
that a dense system was created in the collision whose extreme condition has 
significantly affected particle production mechanisms. Furthermore, various 
observations signal collective (hydro) dynamical behaviour in the collision re- 
gion which in turn indicates the importance of equilibration processes for the 
understanding of the collision dynamics. Especially, the hadronic momen- 
tum spectra show signs of both radial and azimuthally directed flow, and 
two-particle correlations indicate a strong transverse expansion of the source 
before freeze-out. Despite the rich body of these and other observations, it 
remains however controversial to what extent these data are indicative for 
the creation of a QGP or can also be explained in purely hadronic scenarios. 

To make further progress on this central issue, a more detailed under- 
standing of the space-time geometry and dynamics of the evolving reaction 
zone is required. The systems created in relativistic heavy-ion collisions are 
mesoscopic and shortlived, and the geometrical and dynamical conditions 
of the cauldron play an essential role for the particle production processes. 
For example, the maximal energy density attained in the collision, the time- 
dependence of its decrease, and the momenta of the produced particles rel- 
ative to the collectively expanding hadronic system will affect the observed 
particle ratios. Two-particle correlations provide the only known way to ob- 
tain directly information about the space-time structure of the source from 
the measured particle momenta. The size and shape of the reaction zone 
and the emission duration become thus accessible. In combination with the 
analysis of single particle spectra and yields, it is furthermore possible to 
separate the random and collective contributions to the observed particle 
momenta. This permits to also reconstruct the collective dynamical state 
of the collision at freeze-out. These new pieces of information give powerful 
constraints for dynamical model calculations; they can also be taken as an 
experimental starting point for a dynamical back extrapolation into the hot 
and dense initial stages of the collision. The present work reviews the foun- 
dations of HBT interferometry in particle physics and discusses the technical 
tools for its quantitative application to relativistic heavy-ion collisions. 
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1.1 Historical overview 



HBT intensity interferometry was proposed and developed by the radio as- 
tronomer Robert Hanbury Brown in the fifties, who was joined by Richard 
Twiss for the mathematical analysis of intensity correlations. Their original 
aim was to bypass the major constraint of Michelson amplitude interferom- 
etry at that time: in amplitude interferometry, the resolution at a given 
wavelength is limited by the separation over which amplitudes can be com- 
pared. Hanbury Brown started from the observation that "if the radiation 
received at two places is mutually coherent, then the fluctuation in the in- 
tensity of the signals received at those two places is also correlated" [ 
More explicitly, amplitude interferometry measures the square of the sum of 
the two amplitudes Ai and A2 falling on two detectors 1 and 2: 

\A, + = + + {Al A2 + A, a;) . (1.1) 

The last term, the 'fringe visibility' V, is the part of the signal which is sen- 
sitive to the separation between the emission points. Averaged over random 
variations, its square is given by the product of the intensities landing on the 
two detectors [ ^ , 

{V^) = 2 (lAiHAaH + {AfAl) + {AjAf) 2{h h) ■ (1.2) 

The last two terms of this expression vary rapidly and average to zero. Ac- 
cording to ( |1.2|) , intensity correlations between different detectors contain 
information about the fringe visibility and hence about the spatial extension 
of the source. To demonstrate the technique, Hanbury Brown and Twiss 
measured in 1950 the diameter of the sun, using two radio telescopes oper- 
ating at 2.4 m wavelength, and determined in 1956 the angular diameters of 
the radio sources Cas A and Cyg A. Furthermore, they measured in a highly 
influential experiment intensity correlations between two beams separated 



from a mercury vapor lamp. They thus demonstrated [ ^ that photons in 
an apparently uncorrelated thermal beam tend to be detected in close-by 
pairs. This photon bunching or HBT-effect, first explained theoretically by 



Purcell [ 134| |, is one of the key experiments of quantum optics [ How- 



ever, with the advent of modern techniques which allow to compare radio 
amplitudes of separated radio telescopes, Michelson interferometry has again 
completely replaced intensity interferometry in astronomy. 
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In particle physics, the HBT-effect was independently discovered by G. 
Goldhaber, S. Goldhaber, W.Y. Lee and A. Pais [ In 1960, they stud- 
ied at the Bevatron the angular correlations between identical pions in pp- 
annihilations. Their observation (the "GGLP-effect"), an enhancement of 
pion pairs at small relative momenta, was explained in terms of the finite 
spatial extension of the decaying pp-system and the finite quantum mechan- 
ical localization of the decay pions [ ^]. In the sequel of this work, it was 
gradually realized that the correlations of identical particles emitted by highly 
excited nuclei are sensitive not only to the geometry of the system, but also 
to its lifetime [ ^ |15(]|| . This point has become increasingly more impor- 
tant, and it was supplemented by the later insight that the pair momentum 
dependence of the correlations measured for relativistic heavy-ion collisions 
contains information about the collision dynamics [ |125|| . The origins of the 



wide field of applications to relativistic heavy-ion collisions can be dated 
back to the works of Shuryak 



150|| , Cocconi 



Grishin, Kopylov and 
Podgoretskii [ ^ , and to the seminal paper of Gyulassy, Kauffmann 

and Wilson [ Important contributions in the eighties include a more 
detailed analysis of the role of final state interactions [ |9^, ^ |126| , |3^ , ^], 
the development of a parametrization [ 123] taking into account the longi- 

and the 



125, 126 



tudinal expansion of the system created in the collision 
first implementation of the HBT-effect in prescriptions for event generator 
studies [ p.85|| . Also, the effect was seen in and analyzed for high energy colli- 
sions (see the recent review by Lorstad [ |103|| ). In addition, there is a wealth 
of experimental data and theoretical work on correlations between protons 
and heavier fragments (pp, pd, p^He) in lower energy (< 1 GeV) nuclear 
collisions, which are summarized in the review article of Boal, Gelbke and 
Jennings 



32 



With the advent of relativistic heavy-ion beams at CERN and Brookha- 
ven, many of these concepts had to be refined and extended to the rapidly 
expanding particle emitting systems created in heavy-ion collisions. The 
relativistic collision dynamics plays an important role in the derivation of 
the HBT two-particle correlator and of its modern parametrizations. It is 
adequately reflected in recent model discussions of the particle phase-space 
density from which the two-particle correlator is calculated. Several smaller 
reviews [ |103 , 



201, |77|, [T32|, |78 



T70 | 



!2| as well selected reprint volume 
exist by now. The present work aims at a unified presentation of the un- 
derlying concepts and calculational techniques, and of the phenomenological 
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applications of HBT interferometry to the rapidly expanding sources created 
in these relativistic heavy-ion collisions. It does not provide a comprehensive 
review of the experimental data, for which we refer to the overview given in [ 
|8l. 



1.2 Outline 

We start by discussing the relation between the single-particle Wigner phase- 
space density S{x, K) of the particle emitting source, the triple-differential 
one-particle spectrum EpdN/d^p and the two-particle correlation function 
C(q, K) for pairs of identical bosons: 

= fd*xS(x,p), (1.3) 



1 + \{K) exp 



:i.5) 



The approximations are discussed in the main text; the notation used here 
and throughout this review is compiled at the end of this introduction. The 
main aim of particle interferometric methods is to extract as much infor- 
mation as possible about the emission function S{x,K), which characterizes 
the particle emitting source created in the heavy-ion collision. We discuss 
how the above expressions are modified to include final state interactions 
and multiparticle symmetrization effects and how they apply to numeri- 
cal event simulations of relativistic heavy-ion collisions. Contact between 
theory and experiment is made with the help of Gaussian parametrizations 
( |1.5D of the correlator which we review in chapter ^ We discuss the Carte- 
sian Pratt-Bertsch parametrization as well as the Yano-Koonin-Podgoretskii 
(YKP) parametrization where the latter is particularly adapted to the de- 
scription of systems with strong longitudinal expansion. We then turn to 
estimates of the pion phase space density based on such Gaussian fits. Our 
main focus is on the space-time interpretation of the HBT radius parameters 
R^j{K) which we establish in terms of space-time variances of the Wigner 
phase-space density S{x,K). Particle emission duration, average particle 
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emission time, transverse and longitudinal extension of the source as well 
as position-momentum correlations in the source due to dynamical flow pat- 
terns are seen to be typical source characteristics to which identical particle 
correlations are sensitive. While most of our discussion is carried out for 
central collisions, we also review how this framework can be extended to col- 
lisions at finite impact parameter where the HBT radius parameters depend 
on the azimuthal angle of the emitted particles with respect to the reaction 
plane. Furthermore, we discuss more advanced techniques which do not rely 
on a Gaussian parametrization of the correlation function but require better 
statistics of the experimental data. This concludes our review of existing 
analysis tools. 

Chapter ^ is devoted to applications of the presented framework within 
concrete model studies. We introduce a simple but flexible class of models for 
particle emission in relativistic heavy-ion collisions. These are motivated by 
hydrodynamical and thermodynamical considerations and allow to illustrate 
the main techniques discussed before. Different analytical and numerical 
calculation schemes for the HBT radius parameters are contrasted, and we 
explain which geometrical and dynamical model features are reflected by 
which observables. Then we discuss how resonance decay contributions to 
pion spectra modify these calculations, and we compare the results of this 
model with various other model studies, focussing on the qualitative and 
quantitative differences. All these results are flnally combined into an anal- 
ysis strategy for the reconstruction of the particle emitting source from the 
measured one- and two-particle spectra. The method is illustrated on Pb-|-Pb 
data taken by the NA49 Collaboration at the CERN SPS. 



1.3 Notation and conventions 

We use natural units h = c = ks = I- Unless explicitly stated otherwise, 
pairs or sets of N particles are meant to be pairs or sets of identical spinless 
bosons. In particular, we think of like-sign pions or kaons, the most abun- 
dant mesons in heavy-ion collisions. Most of our discussion carries over to 



fermionic particles by replacing the + signs in ( p..4|) and (|1.5| ) by — signs and 
changing from symmetrized to anti-symmetrized A^-particle states whereever 
they appear in derivations. In what follows, we do not mention the fermionic 
case explicitly. 
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Most of our notation is introduced during the discussion. Variables in 
bold face denote 3-vectors. For simpler reference, we list here some of the 
variables used most frequently. 



Pi = (Ei.pi) detected final state particle momenta, on-shell 
m±_=\Jm? + single particle transverse mass 

(f) azimuthal angle of 
y— \ In -^-^^ (reman y) single particle rapidity 

y {italic y) coordinate in configuration space 
Pj simulated particle momenta, e.g. from Monte 
Carlo simulations 
particle momentum operator 
fi simulated particle positions 
ti simulated particle emission times 
K—\{p\ +P2) average pair momentum, off-shell 
M^—^m? + K\ transverse mass associated with K 
$ azimuthal angle of 

rapidity associated with K 
q = {pi — P2) relative pair momentum, off-shell 
(3 — K/K^ velocity of particle pair (approximately) 

C{q,K) two-particle correlation function, also denoted 
by C{pi,p2) 

S{x, K) single-particle Wigncr density, emission function 
p(a;, p) classical phase-space density 
Viip) covariant one-particle spectrum 
7^2(^1,^2) covariant two-particle spectrum 
M normalization of the correlator 
TV number operator 

(7 quantum mechanical wave packet width 
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Chapter 2 



Particle correlations from 
phase-space distributions 



There are numerous derivations of identical two-particle correlations from 
a given boson emitting source. An (over) simplified argument starts from 
the observation that, after weighting the emission points of a two-particle 
Bose-Einstein symmetrized plane wave \E'i2(a;i, a;2, Pi, P2) = i (e*Pi''^i"'"*P2-a:2 
_j_ ^tpi-x2+ip2-xi^ ^ normalized spatial distribution of emission points p{x), 
the two-particle correlator C{q) is given by the Fourier transform of the 
spatial distribution: 



Extracting the spatial information p(x) from the measured momentum spec- 
tra is then a Fourier inversion problem. The solution is unique if we assume 
p{x) to be real and positive. 

Equation ( |2.1|) remains, however, unsatisfactory since it does not allow 
for a possible time-dependence of the emitter and cannot be easily extended 
to sources with position-momentum correlations. Both properties are indis- 
pensable for an analysis of the boson emitting sources created in heavy-ion 
collisions. A sound starting point is provided by the Lorentz invariant one- 
and two-particle distributions for each particle species 
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2 



(2.1) 




(2.2) 
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dN 

^2(pi,P2) = El E2 ^3^^ ^3^^ = El E2 (a+ a+ ap^Qpi) • (2-3) 

These distributions involve expectation values (...) which can be specified 
in terms of a density operator characterizing the collision process. In most 
applications, (...) involves an average over an ensemble of events. The two- 
particle correlation function of identical particles is defined, up to a propor- 
tionality factor A/", as the ratio of the one- and two-particle spectra: 

Cto.P.)-Af „^^'^-p7> . (2.4) 

In section |2.1| , we discuss its normalization as well as the experimentally 
used method of "normalization by mixed pairs". Sections ^]2| and |2.3| then 



deal with two different derivations of the basic relation (|1.4| ) between the 
two-particle correlation function and the Wigner phase-space density. Final 
state interactions and multiparticle symmetrization effects are discussed sub- 
sequently in sections and We conclude this chapter by discussing 
the implementation of this formalism into event generators. 



2.1 Normalization 

The normalization M of the two-particle correlator ( |2.4| ) can be specified 



by relating the particle spectra to inclusive differential cross sections, or by 
requiring a particular behaviour for the correlator ( |2.4| ) at large relative pair 
momentum q. Pair mixing algorithms used for the analysis of experimental 
data approximate these normalizations. 

2.1.1 Differential and total one- and two-particle cross 
sections 



The one- and two-particle spectra ( p.2|/|2.3|) are given in terms of the one- 



and two-particle inclusive differential cross sections as 



VAp) ^ Eip^, (2.5) 
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They are normalized by 



^ V,{p) = (iV), (2.7) 
/^^^^(Pi'P^) = (^(^-1))' (2-8) 

J Il/i 11/2 

where N = J (i^/caj^ak is the number operator. Two natural choices for the 
normalization J\f in (|2.4|) arise [ |113| , |192|, |176| , |191|] by either taking 



directly the ratio of the measured spectra (p.7|) and ( p.8|) , which results in 

Af=l, (2.9) 



or by first normalizing the numerator and denominator of (|2.4| ) separately to 
unity, which gives 

^f=-J^ . (2.10) 

(iV(iV-l)) ^ ^ 

Since in either case Af is momentum- independent, it does not affect the 
space-time interpretation of the correlation function, with which we will focus 
primarily. 

For Af = I, the correlator equals 1 whenever V2{pi, P2)='Pi{pi)Vi{p2) ■ 
Neglecting kinematical constraints resulting from finite event multiplicities, 
one can often assume this factorization property to be valid for large relative 
momenta q. Since at small values of q the correlation function is larger than 
unity, this generally implies {N{N — 1)) > (iV^), i.e., larger than Poissonian 
multiplicity fiuctuations. This is a natural consequence of Bose-Einstein 
correlations. 

The second choice ( |2.10| ) permits to view the correlator as a factor which 
relates the two-particle differential cross section da^^ / d?pid?p2 of the real 
world (where Bose-Einstein symmetrization exists) to an idealized world in 
which Bose-Einstein final state correlations are absent, 

da^^/d'p, d'p2 = C{q, K) da^^/d^p, d'p2 , (2.11) 

without changing the event multiplicities. Such an idealized world is a nat- 
ural concept in event generator studies which simulate essentially the two- 
particle cross sections da^^/d^pi d^p2. They are typically tuned to reproduce 
the measured multiplicity distributions and thereby account heuristically for 
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the effects of Bose-Einstein statistics on particle production; tlie quantum 
statistical symmetrization of the final state, however, is not a part of the 
code. With the normalization ( [^.1(J| ), the factor C{q,K) in ( |2.11| ) preserves 
the total cross sections, cr^^'^°^ = (t5P'*°*. 



2.1.2 Experimental construction of the correlator 

From the data of relativistic heavy-ion experiments, the two-particle corre- 
lator is usually constructed as a quotient of samples of so-called actual pairs 
and 'mixed' pairs or reference pairs. 

One starts by selecting events from the primary data set. Actual pairs are 
pairs of particles that belong to the same event. Reference pair partners are 
picked randomly from different events within the set of events that yielded 
the actual pairs. The correlation function is then constructed by taking 
the ratio, bin by bin, of the distribution Da of these actual pairs with the 
distribution Dr of the reference pairs [ |183| , |184[ , 



^ , . . number of actual pairs in bin (Aq AK) . ^ 

^^(^«-^^) = uumber of actual pairs m sample ' P'^^' 

^ , , » , number of reference pairs in bin (Ao, AK) , 

""(^"■^'f) = number of reference pairs in Lpla 

^ mM- 

The number of reference pairs for each actual pair, the so-called mixing fac- 
tor, is typically between 10 and 50. It has to be chosen sufficiently large to 
ensure a statistically independent reference pair sample while for numerical 
implementations it is of course favourable to keep the size of this sample as 
small as possible. The two-particle correlator constructed in ( |2.14| ) coincides 



with the theoretical definition (|2.4| ) only if the reference pair distribution 
coincides with an appropriately normalized product of one-particle spectra. 
Since both actual and reference pair distributions are normalized to the cor- 
responding total particle multiplicity, this normalization of ( ^.14| ) coincides 
with the normalization M = {N^ /{N{N — 1)). A different construction of 
the correlator from experimental data has been proposed by Miskowiec and 
Voloshin [ |113|| (see also [ |191|| ). Their proposal amounts to a modification of 



the number of pairs in the sample by which (|2.12| ) and (|2.13| ) is normalized 
and coincides with A/" = 1. 
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In general, the reference sample contains residual correlations, which are 
not of physical origin but stem typically from the restricted acceptance of 
experiments. For these residual correlations, corrections can be employed [ 

ml, mi- 



2.2 Classical current parametrization 

How does one calculate the momentum correlations for identical pions pro- 
duced in a heavy-ion collision? The pion production in a nuclear collision is 
described by the field equations for the pion field (f){x) [ |68[|, 

(n + m^^^ix) = J{x) . (2.15) 

This equation is obviously intractable since the nuclear current operator J{x) 
couples back to the pion field and is not explicitly known. The classical 
current parametrization [ |6^ approximates the nuclear current J{x) by a 
classical commuting space-time function J{x). The underlying picture is that 
at freeze-out, when the pions stop interacting, the emitting source is assumed 
not to be affected by the emission of a single pion. This approximation can 
be justified for high event multiplicities [ For a classical source J{x), 
the final pion state is then a coherent state | J) which is an eigenstate of the 
annihilation operator 

ap\J) = zJ{p)\J). (2.16) 

The Fourier transformed classical currents J are on-shell. Using ( p.l6|) , the 
one- and two-particle spectra ( |2.2| ) and ( |2.3| ) are then readily calculated. 
Usually, the classical current is taken to be a superposition of independent 
elementary source functions Jq: 

J{p) = J2e"^'e'P-^'Mp-p^). (2.17) 

If the phases 0j are random, then this ansatz characterizes uncorrelated 
"chaotic" particle emission, and the intercept A in ( |1.5| ) equals one. In more 
general settings, where the phases (pi are not random, the intercept drops 
below unity. One distinguishes accordingly between a formalism for chaotic 
and partially chaotic sources. 
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2.2.1 Chaotic sources 



Chaotic sources are given by a superposition ( p.l7| ) of elementary sources 



Jo, centered around phase-space points Xj, Pi with random phases 0j. The 
corresponding ensemble average specifying the particle spectra ( [^.2[ / [^.3[ ) is [ 



(O) = Y: PnU / d%d'p,pix„p,) / ^{J\0\J), (2.18) 

N=0 4=1"' "'0 

where Pat is the probability distribution for the number N of sources, J2n=o 
Pn = 1 and the normalized probability p{xi,pi) describes the distribution of 
the elementary sources in phase space. A direct consequence of the ensemble 
average ( 2.18| ) is the factorization of the two-particle distribution into two- 



point functions [ 

V2{PUP2) = ^^^f ~2^^^^ (^l(Pl)^l(P2) + |^j(Pl,P2)P) , (2.19) 



SApi,P2) = ,jE,E,{d;^dp,) = ^E,E,{J*{p,)J{p2)), (2.20) 

where (iV(iV - 1))^ = E?^=o ^iv A^(A^ - 1), {N)p = En=oPnN and the 
expectation value in ( p.2q ) is evaluated according to the prescription ( p.l8| ). 
Here, the integer denotes the number of sources, not the number of emitted 
pions. For a Poissonian source multiplicity distribution the prefactor in ( [^.19| ) 
equals unity. In the derivation of ( p.l9| ) a term is dropped in which both final 
state particles come from the same source. This term vanishes in the large 
N limit 



The factorization in (|2.19| ) follows from the commutation relations, once 
independent particle emission and the absence of final state interactions is 
assumed. Due to its generality, it is sometimes referred to as "generalized 
Wick theorem". The emission function S{x,K) which enters the basic rela- 



tion (|1.4|) can then be identified with the Fourier transform of the covariant 
quantity Sj{pi,p2) [ ^ 0. The latter is given by the Wigner transform of 
the density matrix associated with the classical currents, 

2{2Tif 

for which the following folding relation can be derived 

Sj{x,K) = {N)pJd^zd^qp{z,q)So{x-z,K-q), (2.22) 



Sjix,K)= I T^e-^^-^(J*(x + iy)J(x-|y)), (2.21) 
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So{x,K) = j -^^e-^^-yj;{x + \y).Ux-\y). (2.23) 

Here (iV) p is fixed by normalizing the one-particle spectrum to the mean pion 
multiplicity (A^) p = (N) ; the distribution p and elementary source function 
5*0 are normalized to unity. The full emission function is hence given by 
folding the distribution p of the elementary currents Jq in phase-space with 
the Wigner density So{x,K) of the elementary sources. Wigner functions 
are quantum mechanical analogues of classical phase-space distributions [ 
90| . In general they are real but not positive definite, but when integrated 



over X or K they yield the observable particle distributions in coordinate or 
momentum space, respectively. Averaging the quantum mechanical Wigner 
function S{x, K) over phase-space volumes which are large compared to the 
volume {27ch)^ of an elementary phase-space cell, one obtains a smooth and 
positive function which can be interpreted as a classical phase-space density. 

From the particle distributions (|2.19|/pr20D one finds the two-particle cor- 
relator [ [150|, m, |12|, 13, n 



^) = ^ + Sd^xS[x^p,)Sd^yS[y^p,) ' ^^'^^^ 



if the normalization in (|2.4| ) is chosen to cancel the prefactor in ( |2.19|) . 



Adopting instead the normalization prescription ( p.lO|) leads to a normaliza- 
tion of C{q,K) at q — i> oo which is smaller than unity, since ( p.lO| ) is not 
the inverse of the prefactor in ( p.l9|) . 



2.2.2 The smoothness and on-shell approximation 

The smoothness approximation assumes that the emission function has a 
sufficiently smooth momentum dependence such that one can replace 

S{x, K - \q) S{y, K+\q) ^ S{x, K) S{y, K) . (2.25) 

Deviations caused by this approximation are proportional to the curvature of 



the single-particle distribution in logarithmic representation [ ^ and were 
shown to be negligible for typical hadronic emission functions. Using this 
smoothness approximation, the two-particle correlator ( |2.24| ) reduces to the 
expression on the r.h.s. in (|1.^ 
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Eq. ( |1.4| ) which uses the smoothness approximation, forms the basis 
for the interpretation of correlation measurements in terms of space-time 
variances of the source as will be explained in section ^. For the calculation of 
the correlator from a given emission function, the smoothness approximation 
can be released by staring directly from (|2.24| ). In an analysis of measured 



correlation functions in terms of space-time variances of the source, one can 
correct for it systematically [ ^ using information from the single particle 
spectra. 

The emission function S{x, K) depends in principle on the off-shell mo- 
mentum where = + E2). In many applications one uses the 
on-shell approximation 



S{x,K°,K) ^ S{x,Ek,K) , Ek = Vm2 + . (2.26) 

Again, the corrections can be calculated systematically [ ^ but were shown 



to be small for typical model emission functions for pions and heavier hadrons. 
The on-shell approximation (|2.26 ) is instrumental in event generator studies. 



where one aims at associating the emission function S(x, K) with the simu- 



lated on-shell particle phase-space distribution at freeze-out, see section ^[6 
It is also used heavily in analytical model studies, see section ^. 

2.2.3 The mass-shell constraint 



Although the correlator (|2.24|) is obtained as a Fourier transform of the 



emission function S{x,K), this emission function cannot be reconstructed 
uniquely from the momentum correlator (|2.24j) . Note that since the Wigner 



density S{x, K) is always real, the reconstruction of its phase is not the issue. 
The reason is rather the mass-shell constraint 

K-q=pl-pl = ml-ml = Q, (2.27) 

which implies that only three of the four relative momentum components are 
kinematically independent. Hence, the g-dependence of C{q,K) allows to 
test only three of the four independent x-directions of the emission function. 
This introduces an unavoidable model-dependence in the reconstruction of 
S{x,K), which can only be removed by additional information not encoded 
in the two-particle correlations between identical particles. Eq. (|2.27| ) sug- 



gests that this ambiguity may be resolvable by combining correlation data 
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from unlike particles with different mass combinations, if they are emitted 
from the same source. Unlike particles do not exhibit Bose-Einstein correla- 
tions, but are correlated via final state interactions and therefore also contain 
information about the source emission function. In this review, we do not 
discuss unlike particle correlations, although this is presently a very active 
field of research [|9^, ^ |100| , |101| , |168| , [156| , |114|| . It is still an open question to 



what extent a combined analysis of like and unlike particles allows to bypass 
the mass-shell constraint ( p.27| ). 

The consequences of the mass-shell constraint are discussed extensively 
in sections |] and ||. 

2.2.4 The relative distance distribution 

For several applications of two-particle interferometry it is useful to reformu- 
late the correlator ( |1.4| ) in terms of the so-called normalized relative distance 
distribution 

d{x,K) = j d^X s{X + %K)s{X -l,K) , (2.28) 

constructed from the normalized emission function s{x, K). Note that d{x, K) 
= d{—x, K) is an even function of x. This allows to rewrite the correlator 
(0)as 

C(qr, K)-l = J d^x cos(g ■ x) d{x, K) , (2.30) 

where the smoothness and on-shell approximations were used. With the 
mass-shell constraint in the form = (3 ■ q this can be further rewritten in 
terms of the "relative source function" ^^(r): 

C{q,K)-l = J d^x cos{q- x) J dt d{x + f3t,t; K) 

= J d^x cos(q ■ x) SKix) . (2.31) 

In the rest frame of the particle pair where /3 = 0, the relative source function 
Sk{x) is a simple integral over the time argument of the relative distance 
distribution d{x,t;K). In this particular frame the time structure of the 
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source is completely integrated out. This illustrates in the most direct way 
the basic limitations of any attempt to reconstruct the space-time structure 
of the source from the correlation function. 



2.2.5 Partially coherent sources 



It is well-known from quantum optics [ ^ that, in spite of Bose- Einstein 
statistics, the HBT-effect does not exist for particles emitted with phase 
coherence, but only for chaotic sources. This is why in (|2.17|) a chaotic 
superposition of independent elementary source functions Jq{x) was adopted. 
The question of possible phase coherence in pion emission from high energy 
collisions was raised by Fowler and Weiner in the seventies [|^, so far 

the dynamical origin of such phase coherence effects has however remained 
speculative. Their consequences for HBT interferometry can be studied by 
adding a coherent component to the classical current discussed above, 

J{p) = JcoM + JcUp) ■ (2.32) 



An analysis similar to the one presented in section p.2.1| shows that as the 
number ricohip) of coherently emitted particles increases, the strength of the 
correlation is reduced \ |68ll : 



X{K) = 1-D\K), (2.33) 
ricohiK) + ncha(-f^) ' 



D(K) = ^coh(ii:) . . 



For this reason the intercept parameter is often referred to as the 

coherence parameter. In practice various other effects (e.g. particle misiden- 
tification, resonance decay contributions, final state Coulomb interactions) 
can decrease the measured intercept parameter significantly. Although ex- 
perimentally it is always found smaller than unity, < A(i^) < 1, this can 
thus not be directly attributed to a coherent field component. For a detailed 
account of the search for coherent particle emission in high energy physics 



we refer to the reprint collection [ |170|| . Recent work [ |80[ shows that the 



strength D(K) of the coherent component can be determined independent 
of resonance decay contributions and contaminations from misidentified par- 
ticles if two- and three-pion correlations are compared, see section |4.3| . A 
coherent component would also affect the size of the HBT radius parameters 
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and their momentum dependence [ |169| , |Tl|, 0, |80[ . The ansatz ( p.32| ) 
is only one possibihty to describe partially coherent emission. Alternatively, 
one may e.g. choose a distribution of the phases 0j in ( |2.17D which is not com- 
pletely random, thereby mimicking partial coherence 
of these two approaches still remains to be studied. 



148 1 . The equivalence 



2.3 Gaussian wave packets 



It has been suggested repeatedly [ |TT|, |T20|, [TTOl |T7|, |19g, |5T|, [192 



|176|| that due to the smallness of the source in high energy and relativistic 
heavy-ion physics, particle interferometry should be based on finite size wave 
packets rather than plane waves. This leads to an alternative derivation 
of the basic relations ( |1.3D , ( p.. 41 ) which replaces the classical currents from 
the previous section by the more intuitive notion of quantum mechanical 
wave packets, at the expense of giving up manifest Lorentz covariance in 
intermediate steps of the derivation. One starts from a definition of the 
boson emitting source by a discrete set of phase-space points {ri,ti,Pj) or 
by a continuous distribution p(f"j, ti, pj). These emission points are associated 



with the centers of Gaussian one-particle wave packets /«, [ 110, 174, 192 



(^r.Xp,) U{x,U) = . (2.35) 

The wave packets fi are quantum mechanically best localized states, i.e., 
they saturate the Heisenberg uncertainty relation with Axi = cr/V^ and 
Api = l/\/2a for all three spatial components / = 1,2,3. Here, (Ax/)^ = 
{fi, xffi) —{fi, xifi)'^, xi being the position operator, and analogously for Ap/. 
We consider the free time evolution of these wave packets determined by the 
single particle hamiltonian Hq, 

f.{x,t) = (e-^°(*-*^)/,)(a.,t) 

= ■^y'rf3fc/i(fc)e*('=-^-^'^-(*-*^«. (2.36) 

In momentum space, the free non-relativistic and relativistic time evolutions 
differ only by the choices = fc^/2m and = Vk"^ + rn"^, respectively. 
For the non-relativistic case, the integral (|2.36|) can be done analytically. 
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2.3.1 The pair approximation 



Here, we derive the correlator in the so-called pair approximation in which 
two-particle symmetrized wave functions ^ij{x,y,t) are associated with all 
boson pairs constructed from the set of emission points [ |174| , |176| |: 

^ijix,y,t) = -^[f,ix,t)f,iy,t) + f,iy,t) f,ix,t)) . (2.37) 

The norm of this two-particle state differs from unity by terms proportional to 
the wave packet overlaps (/,, fj), but in the pair approximation this difference 
is neglected, {fi, fj) ~ Sij. In section |2]^ we will release this approximation 
and instead start from properly normalized A^-particle wavef unctions. It is 
then seen that the pair approximation is equivalent to approximating the 
two-particle correlator from an A^-particle symmetrized wavefunction by a 
sum of contributions involving only two-particle terms ^ij. 

The two-particle Wigner phase-space density yVij{x,y,pi,p2,t) associ- 
ated with $jj reads [ |9^ 



Wijix,y,pi,p2,t) = I d'x'd'y'<^i,[x + ^,y + ^,t) e'P'-'^' 



X 



^ t 

2 ' 



(2.38) 



Integrating this Wigner function over the positions x, y, we obtain the pos- 
itive definite probability 'Pij{pi,p2,t) to measure the bosons of the state $jj 
at time t with momenta pi, p2. As long as final state interactions are ig- 
nored, this probability is independent of the detection time. For Gaussian 
wave packets it takes the explicit form (the energy factors ensure that Vij 
transforms covariantly) 



El E2 



1 



J d^x d^yWij{x,y,pi,p2,t) 



i2nY 

^Wi{pi,Pi)Wj{p2,P2) + lwi{p2,P2)Wj{pi,Pi) 
+Wi{pi,P2)Wj{pi,P2) 

X COS ((ri-f-j) • (P1-P2) - {ii-ij){Ei-E2 

,2 



Wi{pi,P2) 

sAK) 



P2) exp 



TT cr^ exp 



-K) 



P2) 



(2.39) 
(2.40) 
(2.41) 
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Here, the integral over Si{K) is normalized to unity, and the two-particle 
probability is normalized such that its momentum integral equals one for 
pairs which are well-separated in phase-space. To relate this formalism to 
the emission function ( p.22|/pl23D of the classical current parametrization, we 
rewrite the two-particle probability (|2.39|) in terms of the Wigner densities 
Si{x, K) of the wave packets. 



2 ) 



1 



TT 



: e 

3 6{t - k) e-^-^-r^fl-'-HK-p^f ^ (2.42) 



-E'pi 2 



J d'^x Si{x,pi) J d% Sj{y,p2) 

d^xsi{x,p2) J d'^y Sj{y,pi) 
+^ I d'x s,ix, K) e^'^-^ j A s,{y, K) e'^^'^ 
+^ld''xs,{x,K)e-"'-^J d'ys,{y,K)e''^-y. (2.43) 

In the pair approximation, the two pion spectrum V2{pi,P2) for an event 
with N pions emitted from phase-space points (rj,tj,Pj) is a sum over the 
probabilities Vij of all jN{N — 1) pairs («, j). The corresponding expression 
for a continuous distribution p{ri,ii,Pj) of wave packet centers is obtained 
by an integral over (|2.43|) . Defining 

Vp, = d% d^n dU p(r-i, p,) , jvpi = l, (2.44) 

we find 

^2(Pl,P2) = j VpiVpjV^j{pi,P2) 



J d'^x S^p{x,pi) J dS S„p{y,p2) + 



2 

d^x S^p{x, K) e''^-'' .(2.45) 



The index on S^^p indicates that this emission function is constructed from 
a superposition of wavepackets while the emission function Sj in ( p.22| ) 
was generated from the classical source currents. Similar to ( ^.22| / P^ ), 
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Sv/p{x, K) is given by a folding relation between the classical distribution p 
of wave packet centers and the elementary source Wigner function so(a;, K), 

S^p{x,K) = J d%£fidtip{fi,ti,pi) so{x-ri,t-ti,K-Pi) , {2A6) 

soix, K) = ^ exp f-^ - a^K^\ 5{t) . (2.47) 



The normalization of (|2.46|) is consistent with the interpretation of the inte- 
gral (|1.3| ) as the one-particle spectrum. 

To determine the normalization N{pi,p2) of the two-particle correlation 
function 

we proceed in analogy to the experimental practice of "normalization by 
mixed pairs": An uncorrelated (mixed) pair is described by an unsym- 
metrized product state 

<l>^^^°^^{x,y,t) = Mx,t)f)iy,t), (2.49) 

for which the two particle Wigner phase-space density and the correspond- 
ing detection probability 'Pfj^^°^^{pi,p2) can be calculated [ |174|| according 
to (|38D-(|^. Taking both distinguishable states and ^]f^°'' into 

account and averaging over the distribution p of wave packet centers, the 
normalization N{pi,p2) coincides with the first two terms in ( p.45| ), 

N{pi,P2) = J d'^xS^p{x,pi) J d'^yS^p{y,p2) . (2.50) 

The two-particle correlator then coincides with the basic relation ( |2.24| ) after 
identifying S^p = S. We note already here that starting from a discrete 
finite set of N emission points {{Pi,ri,ti)}i^[i^N], rather than averaging over 
a smooth distribution p, the expression for the two-particle correlator ( p.24| ) 
receives finite multiplicity corrections [ |174 |. These will be discussed in sec- 
tion |2^. 
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2.3.2 An example: the Zajc model 

We illustrate the consequences of the above folding relation ( p.46| / |2lT7| ) with 
a simple model emission function first proposed by Zajc 



S{r,t,p) 

Ms 



Ms exp 



Er, 




(2.51) 



^(27ri?,Po)3' 

This emission function is normalized to a total event multiplicity A^. The 
parameter s smoothly interpolates between completely uncorrelated {s —>■ 0) 
and completely position-momentum correlated {s —>■ 1) sources. In the limit 
s — > 0, this emission function can be considered as a quantum mechanically 
allowed Wigner function as long as RqPq > h/2. In the opposite limit, 

X p 

the position-momentum correlation is perfect, and the phase-space local- 
ization described by the model is no longer consistent with the Heisenberg 
uncertainty relation. Inserting the model emission function ( |2.51| ) into the 
general expression ( p.24| ) for the two-particle correlator one finds [ |186|| 



lim S{x, p) 

s— >1 



5(3) 



m 



(2.52) 



C{q,K) = l + exp(-P|BT9' 

1 



^2 



(2.53) 

- (2W j • 

For sufficiently large s this leads to an unphysical rise of the correlation 
function with increasing q^. One can argue [ |19CI| , |6T| that the sign change in 
( |2.54|) is directly related to the violation of the uncertainty relation by the 
emission function ( p.51| ). 

If one does not interpret ( |2.51 ) directly as the emission function S{x,p), 
but as a classical phase-space distribution p{f,p) of Gaussian wave packets 
centers, then the correlator is readily calculated via ( p.46| ) [ |6T| : 

C{q,K) = l+exp(-p2j^T,a9'' 

1 \ 



R2 
-'^HBT.o- 



1 



P^ = 



Ri + 



a 



{2RPy 



Pn + 



2a2 



(2.55) 
(2.56) 

(2.57) 
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Now 2RP > 1 independent of the value of a, and the radius parameter is 
always positive. Even if the classical distribution p{f,p) is sharply localized 
in phase-space, its folding with minimum-uncertainty wave packets leads to 
a quantum mechanically allowed emission function S{x,p) and a correlator 
with a realistic fall-off in q. 



2.3.3 Spatial localization of wave packets 

Both the two-particle correlator and the one-particle spectrum calculated 
from ( ^.46| ) depend on the initial spatial localization a which is a free para- 
meter. One easily sees that both limits a ^ and a —>■ oo lead to unrealistic 
physical situations: 

In the limit a 0, the wave-packets is sharply localized in coordinate 
space, and the momenta drop out of all physical observables. The one- 
particle spectrum Ep dN/ cPp comes out momentum- independent irrespective 
of the range of the wave packet momenta pi. The momentum correlations 
read [ [l7|, |T9g 



liniC7(pi,p2) = 1 + 



j Vpi Vpj cos [(fi-fj) ■ (P1-P2) - ik-ij) ■ (^1-^2)] • (2.58) 

Due to the cosine term, the dependence of the two-particle correlator on the 
measured relative energy E1—E2 and momentum pi—p2 gives information on 
the initial spatial and temporal relative distances in the source. This is the 
HBT effect. On the other hand ( |2.58| shows that in this limit the correlator 



does not depend on the pair momentum K, since position eigenstates cannot 
carry momentum information. 

In the other limiting case a — > 00, the wave packets are momentum 
eigenstates which contain no information about the emission points rj. In 
this limit, nothing can be said about the spatial extension of the source, since 
the wave packets show an infinite spatial delocalization. A calculation shows 
that also the temporal information is lost in this case: 

jim C(pi,p2) = l + 5pi,p.. (2.59) 

Clearly, physical applications of the wave packet formalism require finite 
values of a. For example, one can use the Gaussian ( |2.51| ) with s = to 
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generate the distribution of wave packet centers. Writing Pq = 2mT to 
allow for an intuitive interpretation of its momentum dependence in terms 
of a non-relativistic thermal distribution of temperature T, the one-particle 
spectrum shows again thermal behaviour ~ exp (— p^/2mTeff), but with a 
shifted temperature [ |110|, |174|, 



Teff = T + -1- . (2.60) 
The corresponding HBT radius parameter reads 

The second terms in these equations reflect the contributions from the in- 
trinsic momenta and spatial extension of the wave packets. 

This shows that repairing possible violations of the uncertainty relation 
in a given classical phase space distribution by smearing it with Gaussian 
wave packets of finite size a, one changes both the single-particle momentum 
spectra and two-particle correlations. While large values of cr strongly affect 
the source size and thereby the HBT-radii but have little effect on the slope 
of the single-particle spectrum, the opposite is true for small values of a. 
With both quantities fixed by experiment, one has therefore limited freedom 
in the choice of a. 

Different attempts to give physical meaning to the parameters a can be 
found in the literature. For example, Goldhaber et al. [ ^ argued that 
the HBT-radius measured in pp annihilation at rest can be interpreted in 
terms of the pion Compton wavelength. Baym recently tried to associate 
a with the coherence length for phase coherence in the source [ On 
the other hand, ( p.60| ), ( p. 611 ) show that (at least in Gaussian models) the 



physical observables have a functional dependence on only two independent 
combinations of the three paramters T, Rq and a. In practice, this allows to 
reabsorb the wave packet width in a redefinition of the source parameters [ 



2.4 Multiparticle symmetrization effects 

Multiparticle symmetrization effects are contributions to the spectra of Bose- 
Einstein symmetrized A^-particle states which cannot be written in terms of 
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simpler pairwise ones. In many-particle systems with high phase-space den- 
sity, the single- and two-particle spectra receive non-negligible contributions 
from multiparticle symmetrization effects. This complicates the interpreta- 
tion of the emission function S{x, K) as reconstructed from the data. 

1301 



T84i g and Pratt [ [128 



Based on strategies proposed by Zajc 
there exists by now an extensive literature on these effects consisting of nu 

|, 118711 and analytical [ |2^, |, |53 



merical 



130, 184 



181, 188, 176, 191 



model studies. Multiparticle symmetrization effects have been considered es- 
sentially in two different settings. Either one starts from events which at 
freeze-out have a fixed particle multiplicity [ p.85| , p.76| , p.89|| , encoded e.g. 
in the model assumptions by choosing sets of phase-space points (pj, rj, ii). 
Bose-Einstein correlations in the final state then lead to an enhancement of 
the two-particle correlator at small relative pair momentum, but they do not 
affect the particle multiplicity. A second approach [ |12^ , |^, |192| , 



does not only calculate the HBT enhancement effect of identical particles, 
but aims at accounting for the effects of Bose-Einstein statistics during the 
particle production processes as well. As a result, modifications of the mul- 
tiplicity distribution of event samples are calculated. 

Here, we first review the formalism for fixed event multiplicities, which 
is tailored to calculate the final state HBT effect only. Then we discuss 
shortly how this formalism can be adapted to calculate changes of multiplicity 
distributions. 



2.4.1 The Pratt formalism 



In his original calculation [ |128| , |132|| , Pratt starts from the (unnormalized) 
probability PAr(p) for detecting A^ particles with momenta p = (pi, . . . ,Pn)- 
It is expressed through single particle production amphtudes Ta{x) for par- 
ticles with quantum numbers a and A^-particle symmetrized plane waves 
U{xi, . . . ,Xn',Pi, ■ ■ ■ ,Pn) as follows: 

-pNip) = (2.62) 
J d^Xi - ■ -d^XNTa^^Xi) ■ ■ ■Taj,{xN)U{xi,. . . ,xn]Pi,- ■ ■ ,Pn) 



{ai} 

U{x, 



.xn;pi, 



,Pn) 



N 



I Pj ■ Xgj 



(2.63) 
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The sum runs over all permutations s of particles. The main assumptions 
entering here are (i) the absence of final state interactions which allows the 
plane wave ansatz ( |2.63D , and (ii) the assumption of independent particle 
emission which allows to factorize the A^-particle production amplitude into 
N one-particle production amplitudes Ta{x). It is technically convenient to 
change from these to the corresponding Wigner transformed emission func- 



tion [ m, |T32 



St{x,p) = E / d'ye'^-yT:{x + f ) Ux - f ) . (2.64) 

Calculating V2{pi,P2), one recovers in this formalism up to a normalization 
factor the usual expression ( p.l9| ) for the two-particle spectrum. This is, how- 
ever, not the relevant calculation because it gives only the two-particle spec- 
trum from events with exactly two particles. The aim of Pratt's formalism 
is to compute the one- and two-particle spectra for events with multiplicity 
A^, including all multiparticle symmetrization effects. They are obtained by 
integrating ( |2.62| / pTBB| ) over N — 1 or N — 2 momenta, respectively. We use 
the notation Vn for iV-particle momentum distributions which, in contrast 
to (2^/2^), are normalized to unity. Using the following building blocks [ 
T28| 



GliPl,P2) 
GniPl,P2) 
Cn 



d^xST{x,K) e"*^-^ 



E2 
d^p 



-Gi{pi,k2)Gi{k2,k3 



Gn{p,p) 



(2.65) 

■G,ikrr,P2) (2.66) 

(2.67) 



one obtains the desired spectra by the following algorithm 
m|T7§: 



185, pa, 38, 51 



N 



; Uj(N — m) GmiP-, P) 



(2.68) 



V%(pi,P2) 



j-i 



N 



J) 



J2 [Gi{pi,Pi)Gj^i{p2,P2) + Gi{pi,p2)Gj-i{p2,Pi)] , (2.69) 
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N 



N 



J2u;{N-m)C„ 



(2.70) 



m=l 



While the sum in (|2.63|) runs over A^! terms, this algorithm involves only sums 
over all partitions {n, In)^ of N elements; this reduces the complexity of the 
problem considerably. The algorithm ( |2.68|) -( p.70| ) is sometimes referred to 
ring algebra" [ |192| 



as 



since the building blocks Cm and Gm have a 
very simple diagrammatic representation in terms of closed and open rings [ 
|128| , |137| , |176|| . It is also referred to as Zajc-Pratt algorithm, since Zajc had 
analyzed essential parts of the above combinatorics in [ |185|| . The definition 
of the Cm sometimes differs by a factor m from the one given here, which 



results in appropriately modified combinatorial factors in ( [^.68 )-( p.70| ). 



While the set of equations (|2.68|) - (|2.70|) constitutes a great simplification 



over a direct evaluation of ( |2.62| ), the high- dimensional integrations required 
to determine Cm in ( p.66|) still limit its applications significantly. Numerical 



Monte Carlo techniques have been proposed [ |128| , |130| , |55| to calculate (p.66|) . 
An alternative strategy can be applied to a small class of simple (Gaussian) 
models, where one can control the m-dependence of Cm analytically or via 
simple recursion schemes. Especially for Gaussian emission functions, ( p.66| ) 
allows for simple one-step recursion relations [ |128 , 
Gn+i and G„ which can be solved analytically [ |192|| . 



130, 39, 1881 between 



2.4.2 Multiparticle correlations for wave packets 



There have been several recent attempts to combine the Pratt formalism with 
an explicit parametrization of the source in terms of A^-particle Gaussian 
wave packets [ |192| , |5T| , |176|| . The strategy in these studies is to associate 
with each event {(?"«, ti,Pi)}ie[i,Af] a properly symmetrized A^-particle wave 
function [ [19|, 0, |1^ 



N 



{{fi,ii,Pi)}ie[l,N] — ^ ^Ar(^,t) 



(2.71) 



where the fi are the Gaussian wave packets of ( p.35| ). Note that the normal- 
ization A/* = l/(\E'Ar|\I^Ar) of n depends on the positions {(r,, tj, Pj)}jg[i^7v] 
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of the wave packet centers in phase-space. As we shall see, this prevents a 
straightforward application of the Pratt formalism. 

The normalized probability Vn{P] {z}) for detecting particles with 
momenta p = (pi, . . . ,Pn) in the specific wave packet configuration {z} = 
{(fj, tj, jv] takes the following form [ |176|| 



(2.72) 



S,s'eVn 1 = 1 



where the building blocks Tij{p) are given in terms of the Fourier transforms 
fi of the single-particle wave packets as follows: 



Ti,{p) = D,{p,t)D*ip,t) 



D^{p,t) 



-iEp{t-ti) 



fi{p,ii) . 



(2.73) 
(2.74) 



The time dependence of ( 2.74 ) drops out in jFj^. From ( |2.72| ) the normalized 
one- and two-particle momentum distributions are obtained by integrating 



over the unobserved momenta [ |185|, |17G 



_ 1 N 

■ s,s'eSN 1=2 



Si ' 



'P%{pi,P2;{z}) = A/"*^ E ^s[s^iPl)^s',S2{P2)l[fs[si 



N 



d^pTij{p) 



(2.75) 

(2.76) 
(2.77) 



Similarly, higher order particle spectra V ^ contain m factors J^s'.Si in each 
term. 

The factors H fij occurring in ( |2.75| / |2r76|) reflect the multiparticle sym- 
metrization effects on the one- and two-particle spectra. They involve the 
overlap between pairs of wave packets /j, /j, which for the simple case of 
instantaneous emission, ii = ij, take the simple form 



I/. 



exp 



\Zi 



-Vj + lap. . 
a 



(2.78) 
(2.79) 
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This overlap equals 1 for i = j and decreases like a Gaussian with increasing 
phase-space distance \zi — Zj \ between the wave packet centers. According to 
( |2.79D , this distance depends on the wave packet width a, and in the limiting 
cases a ^ and a ^ oo, the overlap functions reduce to what is known as 
the 



pair approximation: /j, = 5, 



(2.80) 



As we will see, these limits correspond to the case of infinite phase-space 
volume, i.e., vanishing phase-space density of the source. Then all sums over 
iSat in (|2.75|/pT6D are trivial and the two-particle spectrum ( |2.76|) reduces 
to a sum over all particle pairs, involving only two-particle symmetrized 
contributions, thus coinciding 



17611 with the correlator derived in section ETSl 



In order to apply the Zajc-Pratt algorithm ( |2.68| )- (|2.70|) , the distributions 
given above must be averaged over the phase-space positions (f-j, tj, of the 
wave packet centers: 



TV 



P^(pi,...,p„) 



l[Vp,]V^ipi,...,Pn,{^})- (2.81) 



\i=l 



where Vpi is defined in (|2.44|) . This is the analogue of the sum over quantum 



numbers in ( 2.62 ). At this point, one encounters the problem that the 
normalization jVq, of the A^-particle wave packet does not factorize. This 
destroys the factorization property ( p.66| ) of the Pratt algorithm. However, 
an analytical calculation becomes possible if a different averaging procedure 
is used instead of ( |2.81| ): 



N 



[[Vpi — >Vp = — — [[p{ri,ti,Pi) 



2.82) 



i=l 



(^) = / \^]ld''P^d^r^di^p{r^,kP^)\ (*Tv|^iv) • (2.83) 



Ki=l 



This modification ( p.82| ) is equivalent to working with unnormalized N- 
particle wave packet states, as done e.g. in Ref. [ Zimanyi and Csorgo 
[ p.92|| have tried to give this modification a simple physical interpretation 
by noting that the factor (^fTvl^I/Ar) can be interpreted as an enhanced emis- 
sion probability for bosons (described by normalized wave packets) which 
are emitted close to each other in phase-space. According to (|2.82|) , which 
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no longer factorizes, this version of "stimulated emission" leads to specific 
correlations among the emission points {fi,ti,pj), i.e., the particles are not 
emitted independently. 
For the average 

i UM) , the spectra can be obtained from the Zajc- 
Pratt algorithm (g:68D-(^) using the building blocks 

= [ £pG^{p,p). (2.85) 



One only needs to replace in ( |2.65| ) Pratt's definition ( |2.64D for the single- 



particle Wigner density, St{x^K), by the wave packet analogue S^p{x^ K) 
given in ' ' 



2.4.3 Results of model studies 



Explicit numerical [ p.28| , |130| , |188|| and analytical [ |192| , |176|| calculations 



of multiparticle symmetrization effects have so far only been performed for 
Gaussian source models. In this case the m-th order Pratt terms Gm, see 
(|2.66|), can be calculated analytically. Writing them in the form 



Gm{K+^q,K 



|g)=C^(^?rAPexp 



(2.86) 



the coefficients 



and g^ can then be obtained from simple recursion 
relations [ |192|| . Two generic features are observed in all these studies [ 



1. For increasing m the factors g^^^ become larger. This leads to steeper 
local slopes of the one-particle spectrum for small momenta. Multi- 
particle symmetrization effects thus enhance the particle occupation at 
low momentum. 

2. For increasing m the factors g^"^^ decrease. This broadens the width of 
the two-particle correlator, indicating that multiparticle symmetriza- 
tion effects lead to an enhanced probability of finding particles close 
together in configuration space. 
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While these observations are generic, their quantitative aspects are model 
dependent and sensitive to the particle phase-space density. Writing the 
single-particle spectrum ( |2.68D in the form [ |17(j|| 



^7v(P) = H VrnGrniP,P)/Cm, Vm = ''"^ Cm, (2.87) 

m=i ^y^^) 

the weights Vm, which satisfy Ylm=i = 1, can be analyzed in the limit of 
large phase-space volumes. For a Gaussian source with width parameters Rq 
and Pq in coordinate and momentum space, the phase-space density is given 
by Pps = N/{RqPq)^. In the limit of large phase-space volume {RqPo)^ ^ 1 
one finds for fixed (but not necessarily small) multiplicity [ |176|| : 



^ uj{N-m) 



r ps 



(1 + pps 



(2.88) 



Similarly, the two-particle spectrum can be written as [ |176 



N m-1 

^7v(Pl,P2) = "^^^ ^i'rn-t{Pl,P2) , (2.89) 

m=2 i=l 

Gi{pi,pi) Gm-i{P2,P2) 



Hi^m-i[Pl,P2) 



^ Gi{pi,p2) Gm-»(P2,Pi) ^2 90) 

Ci Crn.—i 



A2;ain the weie;hts normalized to unity, Z]m=2 ""m — 1; same 

limit as above their leading behaviour is given by 

m—2 

Um ^ , , = , , ■ (2.91) 

(l+p™l)- (l + p™l) ^ ^ 

We finally remark, that the correlator obtained from (|2.87| ) and ( |2.89|) takes 
the generic form [ |176| , |191|| 

2 

where A = 1, A/" ~ 1 — (i?o-Po)~'^ for large phase-space volumes, and F{jp\,p2) 
approaches 1 and in the limits q ^ and |q| — > cxd, respectively. 
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2.4.4 Bose-Einstein effects and multiplicity distribu- 
tions 



So far we only discussed multiparticle symmetrization effects for events witfi 
fixed multiplicity N. The question to what extent multipion correlations 
are also reflected in the multiplicity distributions was asked already in the 
seventies [ ^7j. Recently, it was revived in the context of Pratt's formalism [ 
|128| , ^ p.88|| with the aim to calculate the effect of Bose-Einstein statis- 
tics on the particle production process. These applications typically start [ 
128| , |192| , |188|| from a multiplicity distribution p^^ in the absence of Bose- 
Einstein statistics, for example a Poisson distribution = ^ exp[— no] with 
average multiplicity no- For this case the probability pn of finding events with 
multiplicity n after having accounted for Bose-Einstein correlations is then 
computed as 



128, 51 



-1 



\k=0 / 



(2.93) 



where uj{k) is given in ( p. 701) . For this particular multiplicity distribution, the 
multiplicity averaged one- and two-particle spectra are given by the simple 
expressions [ |3^, |192| , |191 



Viip) = H{p,p), (2.94) 
r2{pi,P2) = H{pi,pi)H{p2,P2) + H{pi,p2)H{p2,pi), (2.95) 



where 



oo „ 

H{pi,P2) = E GUpi,P2) = / d''xS{x,K)e'''-\ 

m=l 



(2.96) 



With the effective source distribution S{x, K) introduced in the second step, 
the correlator again takes the simple form ( p,.4D . We expect that this source 
distribution coincides with the one defined in ( |2.21|) in the context of the 
covariant classical current formalism. The reasons are that (i) both satisfy 
( |2.24| and (ii) that the coherent states | J) of ( p.l6|) generate a Poissonian 
multiplicity distribution. 

According to the first equation ( |2.96| ), the emission function S{x,K) 
contains all multiparticle symmetrization effects. Expressed in terms of the 



single-particle Wigner density S'„p in ( |2.46| ), it takes a complicated form. 
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Model studies [ |T2|, |T30|, |3|, |T], |T8|, |T9T[ indicate that irrespective of 



the particular multiplicity distribution the general features discussed below 
( |2.86[ ) persist: compared to the input distribution S^p{x, K), the multiparti- 
cle symmetrized emission function S{x, K) is more strongly localized in both 
coordinate and momentum space. For the intercept parameter A one finds 
results which depend on the specific choice for the multiplicity distribution. 
Cases are known where A decreases strongly with increasing phase-space den- 
sity [|2|,[13|, |188l . 

This discussion illustrates that for sources with high phase-space density, 
where multiparticle symmetrization effects cannot be neglected, the inter- 
pretation of the emission function S{x, K) reconstructed from the one- and 
two-particle momentum spectra by analyzing ( |1.3D -( prBD is not straightfor- 
ward. The question how Bose- Einstein effects on the multiplicity distribution 
and on the phase-space distribution can be disentangled is still open. In the 
remainder of this review we therefore concentrate on the reconstruction of 
S{x, K) and do not further consider its possible contamination by multipar- 
ticle effects. 



2.5 Final state interactions 

Momentum correlations between identical particles can originate not only 
from quantum statistics but also from conservation laws and final state in- 
teractions. 

Energy-momentum conservation constrains the momentum distribution 
of produced particles near the kinematical boundaries. In high multiplicity 
heavy-ion collisions its effects on two-particle correlations at low relative mo- 
menta are negligible. Similarly, constraints from the conservation of quantum 
numbers (e.g. charge or isospin) become less important with increasing event 
multiphcity. Strong correlations exist between the decay products of reso- 
nances, but since resonance decays rarely lead to the production of identical 
particle pairs, they do not matter in practice. 

This leaves final state interactions as the most important source of dy- 
namical correlations. For the small relative momenta q < 100 MeV which are 
sampled in the two-particle correlator, effects of the strong interactions are 
negligible for pions. For protons, however, they dominate the two-particle 
correlations. On the other hand, for pions, the long-range Coulomb interac- 
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tions distort significantly the observed momentum correlations, dominating 
over the Bose-Einstein effect for small relative momenta. Here we discuss 
how Coulomb correlations are calculated for a given source function and how 
they can be corrected for in the data. The aim of Coulomb corrections is 
to modify the measured two-particle correlations in such a way that the re- 
sulting correlator contains only Bose-Einstein correlations, while the effects 
of final state interactions have been subtracted. For this, several simplified 
procedures have been used in the literature, which we review in what follows. 



2.5.1 Classical considerations 



The Coulomb interaction between particle pairs accelerates them relative to 
each other, thus depleting (enhancing) the two-particle correlation function 
at small relative momenta for like-sign (unlike-sign) pairs. In a high multi- 
plicity environment this final state interaction can be reduced by screening 
effects until the particle pair has separated sufficiently from the rest of the 
system. Both these effects can be taken into account in a classical toy model 
which neglects the Coulomb interaction between the pairs for separations less 
than ro and includes it for larger separations [ pl| . The initial and the finally 



observed relative momenta and q are then related by 

iq/2r _ {q,l2f ^ ^ 
2 /i 2 yU ro ' 

where // is the reduced mass. For two pions, = 1X1^,12 and a radius ro 



(2.97) 



10 



fm, this results in a shift {q/2f = {qo/2f ± 20MeV^ The modification of 
the two-particle correlator is then given by the Jacobian \d^qQ/(fq\ = q^jq 
and reads [ |69l BTI 



C{q) 



(Pqa 



(Pq 



Co(go) = C'o(go). 



It 



2fie^ 



\ ^ro(g/2)^ 



(2.98) 



where 6*0(^0) denotes the two-particle correlator in the absence of Coulomb 
interactions. When comparing ( ^.98] ) with the data, the radius ro can be used 
to accomodate for the dependence of the Coulomb final state interaction 
on the source size. This toy model reproduces the qualitative features of 
experimental data surprisingly well but fails to account quantitatively for 
the correct g-dependence of the correlator at very small relative momenta [ 
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2.5.2 Coulomb correction for finite sources 



For a quantum mechanical discussion of final state Coulomb interactions, 
we associate to the emitted particle pairs a relative Coulomb wavefunction, 
given analytically by the confluent hypergeometric function F, 



r(l + irj) e 



-IT]] 1; Z. 



z± = |(gr ± qr • r) = igr(l ± cos( 



(2.99) 
(2.100) 



Here, r = \r\, q = \q\, and 6 denotes the angle between these vectors. The 
Sommerfeld parameter rj = a/{vrei/ c) contains the dependence on the particle 
mass m and the electromagnetic coupling strength e; we write 



± 



Att q/2 



m e 



(2.101) 



where the plus (minus) sign is for pairs of unlike-sign (like-sign) particles. 

To illustrate the influence of a finite source size in a simple case, we 
take recourse to the relative source function SK{f) defined in ( |2.31|) . This 
function describes the probability that a particle pair with pair momentum K 
is emitted from the source at initial relative distance r in the pair rest frame. 
For sources without x-i^-correlations and neglecting the time structure of 
the particle emission process in the pair rest frame, the corresponding two- 
particle correlation for non-identical charged particle pairs reduces to [ ^ ^ 



(2.102) 



Corrections to this expression are discussed in section p.5.4| . For a pointlike 
source Sj^ir) = 5(3) (r), the correlator ( I^TTU^ ) is given by the Gamow factor 
G{rj) which denotes the square of the Coulomb wavefunction ^^°J^{r) at 
vanishing pair separation r = 0, 



r(l + irj) e 



27r?7 



(2.103) 



For a Gaussian ansatz SKif) oc exp [— r^/4i?^], the dependence of the Cou- 
lomb correlations on the size R of the source is then determined via (|2.102| ), 
see Fig. |2.1|. If the particles are emitted with finite separation r, their 
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Coulomb interaction is weaker and the Gamow factor overestimates the fi- 



nal state interaction significantly, see Fig. |2.1| . The source size thus enters 
estimates of the Coulomb correction in a crucial way, and its selfconsistent 
inclusion in the correction procedure can lead to significantly modified source 
size estimates [ |^ . 




0.05 0.1 0.15 0.2 0.25 




0.05 0.1 0.15 0.2 0.25 
Qinv (GeV) 



Figure 2.1: Unlike-sign correlations for < K± < lOOMeV/c and 4.0 <Y< 
5.0. a) Dashed line: Gamow function, corresponding to vanishing source 
radius. Solid line: fit of the NA49 Pb+Pb data to the function ( glQil) . b) 
The same data compared to calculations based on ( p.l02| ) with a spherically 
symmetric Gaussian source, with R = 0.5 fm, 2 fm, 4.6 fm, and 8 fm (from 
top to bottom). The best fit is obtained for R = 4.6 fm. Figure taken from 



2.5.3 Coulomb correction by unlike sign pairs 

Large acceptance experiments can measure like-sign and unlike-sign particle 
correlations simultaneously. The latter do not show Bose-Einstein enhance- 
ment but depend on final state interactions as well. This opens the possibility 
to correct for the Coulomb correlations in like-sign pairs by using the infor- 
mation contained in unlike-sign pairs. The q-dependence of the unlike-sign 
correlations C^^'^{q, K) is often parametrized by a ginv-dependent simple 
function [ |, ^ 

F(ginv) = 1 + {G{r]+) - 1) e-^'-/«« . (2.104) 
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where q^^^ = — {<fY ^'^d Qo is a fit parameter, q^^^ is the square of 
the spatial relative momentum in the pair rest frame, where g° = 0. For 
small ginv, this function approaches the Gamow factor ( ^.1031 ) for a pointlike 



source, while it includes a phenomenological finite-size correction for large 
relative momentum. More recently, one has started to avoid this intermediate 
step by constructing the corrected correlator C^or^H^; ^) like-sign pairs 
directly using bin by bin the experimental data from the measured like- and 
unlike-sign correlators [ |147| , |13| : 

C[--,\q, K) = C^+-Jiq, K) Cj^-Jiq, K) . (2.105) 

In the absence of Bose-Einstein correlations and for pointlike sources, the 
lefthand side of this equation reduces to unity while the righthand side be- 
comes a product of Gamow factors 

This expression provides an estimate for the accuracy of the correction pro- 
cedure (|2.105| ). It deviates from unity by less than five percent for relative 



momenta q > 8-^ [ 147 |. For pions, only the region g < 10 MeV is affected 
significantly, while for the more massive kaons the whole region g < 25 MeV 
shows an error larger than five percent. Calculating the correction factors in 



( |2.105|) for extended sources, this picture does not change since the main dif- 



ference between like-sign and unlike-sign correlations is due to the different 
Gamow factors, and not to the r-dependent confluent hypergeometric func- 
tion in the relative wavefunctions [ p.54|| . As a consequence, one can obtain 
an improved Coulomb correction for heavier particles by dividing out these 



Gamow factors [ |154|| , 



•-^ corr,improvcdly ~ q^^ )G{rj ) ' (.^-iU^j 

This was shown to work with excellent accuracy for a wide range of source 



parameters [ |154 |. 



2.5.4 General formalism for final state interactions 

We now discuss a general formalism for the discussion of the effects of final 
state interactions, starting from an arbitrary two-particle symmetrized wave 
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function ^'(cci, X2i t) which we expand in terms of plane waves 



(2^ (2^ P2, t) (l)p,,p^ixi, X2, t) , (2.108) 

= (/.K(a;,i)0q/2(r,t) = e-^(^i+^^)*e^^^-="e^«-^ (2.109) 

In the last step, we have changed to center of mass coordinates x = ^{xi+X2) 
and relative coordinates r = {xi — X2). The two-particle state ^' is evolved 
with the interacting hamiltonian H while the plane waves in which we expand 
follow a free time evolution, determined by Hq, Hq, 

Here M = 2m and /j, = m/2 are the pair and reduced mass, respectively. 
The two-particle state ^ determines the two-particle Wigner phase-space 
density and hence the two-particle correlator. The probabihty V^{pi,P2,t) 
of detecting the bosons at time t with momenta pi, p2 is 

'P^{puP2,t) = A%{puP2,t) Am{puP2,t) . (2.111) 

Let us assume that from a time to onwards final state interactions have to be 
taken into account in the description of the time evolution of ^. The time 
evolution of ^*(pi,P2,t) then reads 

A^{puP2,t) = J d'xd'rr^ix,to) [e'(^°+^i-)(*-*°)e-^°(*-*°)0,/2(r,io) 

x^{x + ^,x-^,to). (2.112) 

We are interested in the limit t — > 00 of this expression. To this end, we use 
the M0ller operator 

n+ = lim e^(^o+^rnt)(*-*o) e-iHo{t-to) ^ (2.113) 

t—^0O 
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which determines the solution of the Lippmann-Schwinger equation for the 
corresponding stationary scattering problem, 

(r) = f]+0,/2(r,to). (2.114) 

Irrespective of the form of V{r), once is determined, the two-particle 

detection probability V^{pi,p2,t = oo) is known from (|2.111D . The corre- 



sponding two-particle spectrum is then given by summing over all pair wave 
functions of the event: 

■^2(^1,^2)= J2 Al{pi,p2,t = oo)A^{pi,p2,t = oo). (2.115) 

pairs Vt" 

Coulomb correlations for instantaneous sources 



We now illustrate the use of the two-particle spectrum ( |2.115| ), starting from 



the Gaussian wavefunction introduced in section |2^. The sum Impairs "i- 
( |2.115|) is then a sum over all pairs (i,j) of the set {fi,ti,Pf), or an average 



over some distribution p{ri,ti,Pj). We restrict the calculation to instan- 
taneous emission at time ti = tj = to- For non-identical particles (e.g. 
unlike-sign pions) the two-particle wave function at emission reads then 

'^ij{r,X,to) = ^rcl(r",to) ^pair(aJ,to) 

= (7r(T2)-3/2e-^(*-)'+2i^:-e-3^(-r-)^+§9-r-^ ^2.116) 

and the corresponding amplitude entering the two-particle spectrum ( |2.115| ) 
is 

^*(Pl,P2,t = 00) = (0i^|*pair)($^;/f l^rcl) (2.117) 

t* 



^'^/fl^rei) = /dV$^;/f*(r)vI/rei(r), (2.118) 



(•fxl^pair) = (vra2)-3/V^^(^-^)^e2--(^-^). (2.119) 

For pairs of identical charged particles, the state ( p.ll6|) and the correspond- 
ing amplitude (|2.117|) must be symmetrized properly, adding the missing 
q ^ —q terms and replacing e.g. $q/2 by ± ^-9/2)- Further analyt- 

ical simplifications of the amplitude ( |2.117| ) depend on the functional form 
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assumed for the two-particle state \E'. A study with Gaussian wave packets 
was presented in Ref. [ p.79|| . In the plane wave limit cr — > one recovers 
Pratt's result [ |T5§ 



V^{q,K) oc G{T])[F{-ir];l;iz^)F*{-ir];l;iz_) 
+F{-ir]; l;iz+) F*{-ir]; l;iz+) 
±(;^Qrcosep^_-^. i-iz__)F*{-ir];l;iz+) 
±(,-iircosep^_-^. ^. .^^^ F*{~i7]- 1; iz^)] . (2.120) 

The first two lines are the Born probabilities |$q/2(*")P and |$_q/2(''")P; the 
exchange or interference terms in the last two lines exist only for pairs of 
identical particles. Weighting ( |2.120| ) with the distribution ^^(r), one ob- 



tains the properly symmetrized generalization of ( |2.102|) for pairs of identical 
bosons. 

Coulomb correlations for time-dependent sources 

In general, identical bosons interfering in the final state of a relativistic heavy- 
ion collision are produced at different emission times 7^ tj. This temporal 
structure is neglected in the ansatz for the two-particle wavefunction ( |2.116| ) 
and does not appear in the corresponding result ( |2.12CI| ). A formalism appro- 
priate for the calculation of two-particle spectra from arbitrary space-time 
dependent emission functions S{x, K) was developed in Ref. [ Q (for a rela- 
tivistic approach using the Bethe-Salpeter ansatz see Ref. [ p9|): 



^^^rf^x d'^yS{x,K + p) 



^W,{x-y,p)S{y,K-p), (2.121) 
W,„{x-y,p) = /^e-«-(^-^)x,/2(p+f)x;/2(p-f) .(2.122) 

Here Xq/2 denotes essentially a Fourier transformed relative wavefunction 
times a propagator [ ^, and the function Wq/2 can be interpreted as the 
Wigner density associated with the (symmetrized) distorted wave describ- 
ing final state interactions. For a free time-evolution in the final state, the 
two-particle spectrum ( p. 1211) coincides with the appropriately normalized 



spectrum ( |2.19| ). From the result ( |2.121| ) for the general interacting case. 
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a simplified expression can be obtained by expanding the temporal compo- 
nent of the phase factor in Wq/2 in leading order of the small energy transfer 
caused by the final state interaction [ M 



^2(Pl,P2) 



lP2 



± 



V 



j d^xd\s(x+l,K) S(x-1,k) 

X '^-^/2{y-vy°) ^^/2{y-vy'^) , (2.123) 

(2.124) 



K 



E 



Pi 

Ek 



V2 



P2 

Ek 



The velocities Vi, V2, and v are associated with the observed particle mo- 
menta pi, P2, and their average K. In all three cases the argument of the 
FSI-distorted wave $ can be understood as the distance between the two 
particles in the pair rest frame at the emission time of the second particle. 
Equation ( |2.123| ) is obtained without invoking the smoothness approxima- 



tion. Employing also the latter, both terms in (|2.123|) are associated with the 
same combination of emission functions which can then be written in terms 
of the (unnormalized) relative distance distribution, see also (|2.28|) , 



D{y, K) = j d^x S(x+1,k) S(x- f , K 



(2.125) 



This function denotes the distribution of relative space-time distances y be- 
tween the particles in pairs emitted with momentum K. A particularly simple 
expression due to Koonin is then obtained [ 0, |3^, |[ in the pair rest frame, 
v = Q = K: 



d'y 



dy^Diy,K) 



(2.126) 



where = ':^{^q/2 ± '^-q/2) for identical particle pairs. As explained in 
section |2.2.4| , the last factor in (|2.126|) coincides up to normalization with 
the relative source function ^^(r). With the help of the smoothness ap- 
proximation the two-particle spectrum can thus be expressed by the relative 
source function weighted with the Born probability of the Coulomb relative 
wavefunction, as given before in (|2.102|). 
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We finally mention that first attempts have been made to include in the 
analysis of Coulomb final state effects the role of a central Coulomb charge [ 
IT3, |149| or effects due to high particle multiplicity [ . It is an important 



open question to what extent these effects modify the analysis presented here. 



2.6 Bose-Einstein weights for event genera- 
tors 

Numerical event simulations of heavy-ion collisions provide one important 
method to simulate realistic phase-space distributions. Many such event gen- 
erators exist nowadays. In principle, their output should be a set of observ- 
able momenta Pi with all momentum correlations (and hence the complete 
space-time information) built in. However, none of the existing event gen- 
erators propagates properly symmetrized A^-particle amplitudes from some 
initial condition. As a consequence, the typical event generator output is a 
set of discrete phase-space points (fj,tj,Pj) which one associates with the 
freeze-out positions of the final state particles. This simulated event infor- 
mation {vi, ti,Pj) lacks correlations due to Bose-Einstein symmetrization and 
other types of final state effects. 

We first discuss different schemes used to calculate a posteriori two- 
particle correlation functions for inputs of discrete sets of phase-space points 
{ri,ti,Pj). We then turn to so-called shifting prescriptions which aim at 
producing modified final state momenta with correct particle correlations. 



2.6.1 Calculating C(g, K) from event generator output 

The conceptual problem of determining particle correlations from event gen- 
erators is well-known [ 104| , |l| |105 |: Bose-Einstein correlations arise from 



squaring production amplitudes. They hence require a description of produc- 
tion processes in terms of amplitudes. Numerical event simulations, however, 
are formulated via probabilities. This implies that various quantum effects 
are treated only heuristically, if at all. Especially, event generators do not 
take into account the quantum mechanical symmetrization effects. In this 
sense, the event generator output is the result of an incomplete quantum 
dynamical evolution of the collision. The aim of Bose-Einstein weights is to 
remedy this artefact a posteriori by translating the phase-space information 
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of {fi,ii,pj) into realistic momentum correlations. For a set of A^^ev events of 
multiplicities Nm, this implies formally 

{{{fi,kp,)he[i,N^]] =^ C{q,K) . (2.127) 

The set {{vi, U, Pj)}ie[i,Ar^] denotes the phase-space emission points of the Nm 
like-sign pions generated in the m-th simulated event. The event generator 
simulates thus a classical phase-space distribution 



Pciass(p, r,t) = ^J2T. '^^'H^ - S^'\p - Pd 5{t - k) . (2.128) 

^^ev rn=l i=l 

Prescriptions of the type ( p.l27| ) are not unique: a choice of interpretation is 
involved in calculating two-particle correlations from the event generator out- 
put. Here, we mention two different interpretations of {ri,ii,p^), sometimes 



referred to as "classical" and "quantum" [ ^ 



"Classical" interpretation of the event generator output 



In the "classical" interpretation [ |190| , |178| , |61[] the distribution of phase-space 



points {{(?*«, ti, Pi) }i6[i,Afm]}me[i,Afcv] is interpreted as a discrete approxima- 
tion of the on-shell Wigner phase-space density S{x,p), 

Six,p) = pd^sip,x,t). (2.129) 

The emission function is thus a sum over delta functions. For practical ap- 
plications, it is convenient to replace the delta functions in momentum space 
by rectangular 'bin functions' [ |19CI|| or by properly normalized Gaussians [ 
|178| , pT| of width e — > (we denote both choices by the same symbol Sp^^p) 

_ / 1 : Pi - f < P^,3 <Pl + l U = (n 1 on^ 

'^P^'P ~ \ : else, ^^'^'^^^ 
41 = (J)37^-P(-(P^-P)V^^)- (2.131) 
The one-particle spectrum and two-particle correlator then read [ |190| , |178 



61 



rlN r 1 A^cv A^m 

= / d'-Six,p) = E E41^ (2-132) 



ev rn=l i=l 
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=1 "Pj,pi 



.(e) 

1 "P, >P2 



(2.133) 



The correlator (|2.133|) is the discretized version of the Fourier integrals in 
( |2.24| ). It does not invoke the smoothness approximation, in contrast to the 
popular earlier algorithm developed by Pratt [ |129| , |133|| (which includes final 
state interactions). The subtracted terms in the numerator and denomina- 
tor remove the spurious contributions of pairs constructed from the same 
particles [ |174 . 

In general the result for the correlator at a fixed point (q, K) will depend 
on the bin width e. Finite event statistics puts a lower practical limit on 
e. Tests have shown that accurate results for the correlator require smaller 
values for e (and thus larger event statistics) for more inhomogeneous sources. 
In practice the convergence of the results must be tested numerically 



"Quantum" interpretation of the event generator output 



In the "quantum" interpretation [ 174, 1^1], FTE, 178, Wi] the event generator 
output Pciassir,i,p) is associated with the centers of Gaussian wave packets 
(p.35|) . Neglecting multiparticle symmetrization effects, the corresponding 
Wigner function according to (|2.46| ) reads 

S{x, K) = j £pi Sfi dti pciass(r"i, ii, Pi) so{x-fi, t-ii, K-p^) . (2.134) 

The one- and two-particle spectra are [ |174| , ^ 

dN 



(2.135) 



'ev m=l i=l 



C{q,K) = 1 + 6-^^'"'/'- 



y-A^ev 






y^Afev 


l^miPa) Vm{Ph) - Y.f=l Si{Pa) Si{Ph) 





(2.136) 

Again, the terms subtracted in the numerator and denominator are finite 
multiplicity corrections which become negligible for large particle multiplic- 
ities [ [171 . 
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Discussion 



In both algorithms, the particle spectra are discrete functions of the input 
(fjjtjjpj but they are continuous in the observable momenta pi, p2 and 
hence, no binning is necessary. Each of the sums in ( |2.133D and ( p^.l35| ) 



requires only 0{Nm) manipulations. However, once final state interactions 
are included, the number of numerical operations increases quadratically with 
Nm since the corresponding generalized weights [ do no longer factorize. 
When also accounting for multiparticle symmetrization effects, more than 



0{N^) numerical manipulations are typically required [ |176 



The "classical" and "quantum" algorithms then differ in two points: 

1. There is no analogue for the Gaussian prefactor exp (— cr^ of 
(|2.136 ) in the "classical" algorithm. This is a genuine quantum ef- 



fect stemming from the quantum mechanical localization properties of 
the wave packets. 

2. For the choice a = 1/e, the bin functions 6p^^p are the classical coun- 
terpart of the Gaussian single-particle distributions Si{p). Finite event 
statistics puts a lower practical limit on e, but in the limit e ^ the 
physical momentum spectra are recovered. In contrast, in the "quan- 
tum" algorithm a denotes the finite physical particle localization. In 
this case, the limit cr — > cxd (corresponding to e — > 0) is not physically 
relevant: it amounts according to an emission function with infinite 



spatial extension, yielding limg-^oo C{q, K) = 1 + dq^ [ |174 | 



These algorithms have been shown to avoid certain inconsistencies aris- 
ing from the use of the smoothness approximation for sources with strong 
position- momentum correlations [ |107| , |190|| . Systematic studies indicate that 



violations of the smoothness approximation occur only for emission functions 
S{x, K) which are inconsistent with the uncertainty relation, i.e., which can- 
not be interpreted as Wigner densities. Pratt has shown that typical source 
sizes in heavy-ion collisions are sufficiently large that this problem can be 
neglected [ |133| . 



Extensions of these algorithms to include final state interactions [ H and 



multiparticle correlation effects [ |176|| were proposed but have not yet been 
implemented numerically. 
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2.6.2 Shifting prescriptions 

In the previous subsection we reviewed algorithms which calculate two-parti- 
cle correlation functions from a discrete set of phase-space points. The output 
of the algorithm is a correlator C{q,K) which denotes the probability of 
finding particle pairs with the corresponding momenta; it is not a set of new 
discrete momenta pj with the correct Bose-Einstein correlations included. 
The latter is of interest e.g. for detector simulations which require on an 
event-by-event basis a simulated set of particle tracks to anticipate detector 
performance. Also, it could be used to investigate eventwise fluctuations 
which is not possible with an ensemble averaged correlator. 

The most direct way to achieve this goal would seem to use symmetrized 
amplitudes for the particle creation process. Such a scheme has been devel- 
oped in the context of the Lund string model [ ^, ^ for the hadronization 
of a single string. For more complicated situations there exist so far only 
algorithms which shift after particle creation the generated momenta pj to 
their physically observed values Pj 

Pj =^ Pj- (2.137) 

While the function C{q,K) describes the two-particle correlations only for 
the ensemble average, the set {pj} represents all measurable momenta of a 
simulated single event with realistic Bose-Einstein correlations. 

Such a shifting prescription which employs the full phase-space informa- 
tion of the simulated event was developed by Zajc [ |184| , |185|| . In Ref. [ 



185| | a self-consistent Monte-Carlo algorithm is used for a simple Gaussian 
source to determine the shifts ( p.l37| ) by sampling the momentum-dependent 
iV-particle probability. Although technically feasible, this calculation of N- 
particle symmetrized weights involves an enormous numerical effort. 

Another class of algorithms is used in event generators for high energy 
particle physics. Unlike ( ^.137| ), they explicitly exploit only the momentum- 
space information of the simulated events. Additionally, an ad hoc weight 
function is employed which one may relate to the ensemble-averaged space- 



time structure of the source [ p.04| , |105 



g^^^'^^'"'Q + 5Q. (2.138) 

This shifting procedure involves only particle pairs and is significantly sim- 
pler to implement numerically. By decoupling the position and momentum 
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information one looses, however, possible correlations between the particle 
momenta and their production points. Also, in an individual event this 
shifting prescription is insensitive to the actual separation of the particles in 
space-time. A further problem is that the translation of 5Q from ( p.l38| ) into 
a change of particle momenta is not unique. It changes the invariant mass 
of the particle pair and does not conserve simultaneously both energy and 
momentum. These deficiencies are repaired by a subsequent rescaling of mo- 
menta; according to Ref. [ |105|| the results show in practice little sensitivity 
to details of the implementation. 

A more sophisticated method [ |5^, |5^ attempts to implement the full 
A^-body symmetrization by a cluster algorithm. Again the weights used in 
this algorithm, at least in the present version, encode the space-time structure 
of the source only via a single ensemble-averaged radius parameter. 
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Chapter 3 

Gaussian parametrizations of 
the correlator 



In practice, the two-particle correlation function is usually parametrized by a 
Gaussian in the relative momentum components, see e.g. ( |1.5| ). In this chap- 
ter we discuss different Gaussian parametrizations and establish the relation- 
ship of the corresponding width parameters (HBT radii) with the space-time 
structure of the source. 

This relation is based on a Gaussian approximation to the true space-time 



dependence of the emission function [ 0, |4^, |171| , [79 



S{x,K) = N{K) S{x{K),K)exp 
+6S{x,K). 



(3.1) 



For the present discussion, we neglect the correction term 5S{x,K). We 
discuss in chapter ^ how it can be systematically included. The space-time 
coordinates x^ in (|3.1| ) are defined relative to the "effective source centre" 



x{K) for bosons emitted with momentum K [ |89|, |171| , p3 
x'^(if) = - x^(ir) , x^(K) = (x^)(ir) 



where (. . .) denotes an average with the emission function S'(x, K): 

Jd^xf{x)S{x,K) 



{f){K) 



J d^x S{x, K) 



(3.2) 



(3.3) 
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The choice 

{B-')^,iK) = {x^x,)iK) (3.4) 

ensures that the Gaussian ansatz has the same rms widths in space-time 
as the full emission function. Inserting ( |3.1| ) into the basic relation ( |1.4| ) one 
obtains the simple Gaussian form for the correlator 

C(q, K) = l + exp [-g^g,(£'^F)(K)] . (3.5) 

This involves the smoothness and on-shell approximations discussed in chap- 
ter ^ which permit to write the space-time variances (x^Xy) as functions of 
K only. Note that the correlator depends only on the relative distances x'^ 
with respect to the source center. No information can be obtained about the 
absolute position x{K) of the source center in space-time. 

According to ( p.5|) the two-particle correlator provides access to the rms 
widths of the effective source of particles with momentum K. In general, 
these width parameters do not characterize the total extension of the colli- 
sion region. They rather measure the size of the system through a filter of 
wavelength K. In the language introduced by Sinyukov [ [L53|| , this size is 
the "region of homogeneity" , the region from which particle pairs with mo- 
mentum K are most likely emitted. Space-time variances coincide with total 
source extensions only in the special case that the emission function shows 
no position-momentum correlation and factorizes, S{x,K) = f{x)g{K). 

Relating (|3.5|) to experimental data requires first the elimination of one 
of the four relative momentum components via the mass-shell constraint 
( |2.27|) . Depending on the choice of the three independent components, dif- 
ferent Gaussian parametrizations exist. In what follows, we focus on their 
interpretation in terms of the space-time characteristics of the source. 



3.1 The Cartesian parametrization 



The Cartesian parametrization [ |123| , |125| , p5| , p| is expressed in the out- 
side-longitudinal (osl) coordinate system, defined in Fig. It is based 
on the three Cartesian spatial components Qo (out), Qs (side), qi (long) of 
the relative momentum q. The temporal component is eliminated via the 
mass-shell constraint ( |2.27| ) 

^° -- f3q, f3 = K/K\ (3.6) 



/3 



(/5±,o,a: 



in the os/-system. 



(3.7) 
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This leads to a correlator of the form ( |1.5| ) with X^jj running over j = o, s, I. 
In general this correlator C{q, K) depends not only on K± and Ki, but also 
on the azimuthal orientation $ of the transverse pair momentum | \ . This 
angle, however, does not appear explicitly in the os/-system which is oriented 
for each particle pair differently by the angle $ in the transverse plane. $ 
has to be defined with respect to some pair-independent direction in the 
laboratory system, e.g. relative to the impact parameter b: 

$ = Z(Kx,6). (3.8) 

In the following we discuss both the azimuthally symmetric situation with 
impact parameter 6 = 0, when all physical observables are independent, 
and the parametrization for finite impact parameter collisions. 




Figure 3.1: The os/ coordinate system takes the longitudinal {long) direction 
along the beam axis. In the transverse plane, the out direction is chosen par- 
allel to the transverse component of the pair momentum -fCj,, the remaining 
Cartesian component denotes the side direction. 



3.1.1 Azimuthally symmetric collisions 

For central collisions, 6 = 0, the collision region is azimuthally symmetric, 
and the emission function and correlator are ^-independent. In the osl- 
system, this $-invariance results in a reflection symmetry with respect to 
the side-direction [0: 

^iab(a:; K^, $, Ki) = 5iab(x; K^, $ + 5$, K{) 

<^=^ ^os^a;, y, z, t; Kj_, Ki) = Sos\{x, -y, z, t; Kj_, Ki) . (3.9) 
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We follow common practice in dropping now the subscript osl. Whenever 
azimuthal symmetry is assumed, the emission function S{x, K) is specified 
in the os/-system. 

Due to the y —y reflection symmetry of the emission function, y = 
(y) = and the three space-time variances {x^Xu){K) linear in y vanish. The 
symmetric tensor B^^^K) then has only seven non- vanishing independent 
components. These combine to four non-vanishing HBT-radius parameters 
B^{K) which characterize the Gaussian ansatz ( |1.5|) for the correlator [ 



Rl{K) = {f)iK), (3.10) 

RliK) = {{S:-pJf){K), (3.11) 

R',{K) = {{z - (3itf){K) , (3.12) 

Rl{K) = {{x-f3j){z-f3ii))iK), (3.13) 

RliK) = 0, (3.14) 

RUK) = 0. (3.15) 



Obviously, the more symmetries are satisfied by the emission function, the 
simpler are the expressions obtained for the HBT radius parameters. In 
this context we mention the case of a longitudinally boost-invariant source 
showing Bjorken scaling. Though not strictly satisfied by the finite sources 
created in heavy-ion collisions, this can provide a simple intuitive picture 
of the collision dynamics near mid-rapidity. Longitudinal boost-invariance 
implies a 5 — >■ —z reflection symmetry of the emission function. Thus, in 
addition to the space-time variances linear in y, now also those linear in z 
vanish, and one is left with only 5 non- vanishing independent components of 
B^^{K). In the longitudinally comoving system (LCMS), where Pi = 0, this 
leads to the further simplications 

RUK) = {z'){K) , (3.16) 
RUK)=0, (3.17) 
for longitudinally boost-invariant sources in the LCMS. 

The general relation between the symmetries of the system, the number of 
its independent non-vanishing space-time variances, and the number of non- 
vanishing observable HBT parameters is summarized in the following table: 
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symmetry 




R,,{K) 


none 


10 indep. fcts. 


6 indep. fcts. 






ofir±, 0, F 


azimuthal 


7 indep. fcts. 


4 indep. fcts. 




oiK^, Y 


of i^x, Y 


azimuthal 


5 indep. fcts. 


3 indep. fcts. 


+ long, boostinv. 






in the LCMS 







In all cases, there are more independent space-time variances {x^Xi,){K) 
than experimental observables. This arises from the mass-shell constraint 
( |2.27| ) which leads to a mixing of spatial and temporal variances in the ob- 
servable HBT parameters. One of the most important questions is therefore 
which other properties of the expanding system can be exploited to further 
disentangle spatial and temporal information about the emission function. 

For this, we note that independent of the particular emission function, no 
direction is distinguished in the transverse plane for = 0. The out- and 
s2(ie-components of all observables coincide in this kinematical limit. For the 
space-time variances, this implies 



{zx) 



K^=0 
Ki =0 



if) 

{ix) 



K, =0 



Ki =0 



0. 



(3.18) 
(3.19) 



As long as the limit K_i_ — > of the emission function S{x, K) results in an 
azimuthally symmetric expression (an exception is the class of opaque source 
models discussed in section [5.1.2| ), the above relations between the space-time 
variances at = imply that the HBT radius parameters satisfy 



lim Rl(K) 
lim Rl(K) 



lim RtiK) 
0. 



(3.20) 
(3.21) 



In phenomenological HBT analyses one very often exploits these relations by 
distinguishing between an implicit li'-dependence (due to the X-dependence 
of the space-time variances) and an explicit one (resulting from the mass-shell 
constraint = q ■ K/K^). If the emission function features no position- 
momentum correlation, then all space-time variances are 1^-independent and 
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Eqs. ( |3.18|/prT9D hold independently of K. The difference between R^{K) 



and Rl{K) at non-zero K is then only due to the explicit l^T-dependence 
in (g^TOj/grg), i.e. the term I3\{P). This implies that the explicit K 
dependence dominates if the emission duration is sufficiently large [ |138 



or if the position-momentum correlations in the source are sufficiently weak [ 

Rl{K) - Rl{K) ^ . (3.22) 

In this case, the difference between these two HBT radius parameters gives 
direct access to the average emission duration (t^) of the source and allows 
to partially disentangle the spatial and temporal information contained in 

The gives rise to the following simple interpretation of the HBT radius 
parameters from the Cartesian parametrization: Rg measures the width of 
the emission region in the side direction, and Ro measures the correspond- 
ing width in the out direction plus a contribution from the emission dura- 
tion which can be extracted according to (|3.22| ) under the assumption of a 



weak i^x-dependence of the emission function. The longitudinal radius Ri 
finally describes the longitudinal extension of the region of homogeneity in 
the LCMS where Pi = 0. No easy intuitive interpretation exists for the out- 
longitudinal radius parameter i?^;. It is perhaps best understood in terms of 



the linear correlation coefficient \ 144 



which can be positive or negative but is bounded by \poi{K)\ < 1 due to 
the Cauchy-Schwarz inequality. This coefficient was shown [ p.44|| to be of 
kinematical origin and useful for the interpretation of the longitudinal mo- 
mentum distributions. We shall see in section ^]2| that Roi plays a crucial 
role in the determination of the longitudinal velocity of the emitting source 
volume element. 



3.1.2 Collisions with finite impact parameter 

If the azimuthal symmetry of the particle emitting source is broken, then the 
transverse one-particle spectrum depends on the azimuthal direction of the 
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emitted particles. This can be quantified in terms of tlie liarmonic coefficients 

dN 



dN 



E 



d^p 



E- 

Pt dpt dy d 

E d^N 
27r pt dpt dy 



d xS{x,p) 



l + 2^VniPt,y) cosn{(f) - ipn) 



(3.24) 



n=l 



The size and momentum dependence of the lowest of these harmonic co- 
efficients has been analyzed experimentally at both AGS and CERN SPS 
energies [ [15|, |1^, |1 18| , |124|| . This allows to determine the orientation of 
the reaction plane for semiperipheral collisions event by event with an uncer- 
tainty of less than 30° [ 165 , 124 |. Several attempts to extend this azimuthally 
sensitive analysis to two-particle correlation functions exist [ p.67| , |56| , |175|] . 
The correlation measurements depend on the azimuthal direction $ of the 
pair momentum, see ( p.8|) , and hence allow to provide additional azimuthally 
sensitive information. The corresponding Gaussian radius parameters can be 
written formally in terms of space-time variances which are rotated via P$ 
from the impact parameter fixed to the osl coordinate system 



(/3±, 0, A) • 



(3.25) 



We here differ from the notation adopted in the rest of this review: the 
coordinates x, y and z are here given in the impact-parameter fixed system, 
not the osl one. As for the azimuthally symmetric case, the HBT radius 
parameters show implict and explicit /^-dependences. Their ^-dependence 
thus has two different origins [ |167| , |175| : 

Rl{Kx, $, Y) = (x^) sin^ $ + (f) cos^ $ - (xy) sin 2$ , 
Rl{K^, $, Y) = (x^) cos^ $ + (f) sin^ $ + /3l{P) 

-2/3_L {ix) cos $ - 2/3± (ty) sin $ + {xy) sin 2$ , 
Rl^{K^,^,Y) = {xy) cos2<l> + |sin2$((y2) _ (^2)) 

+(3± {ix) sin $ - {iy) cos $ , 
RUK^,<l>,Y) = {{z-f3iif), 

RliiK^, $, r) = ((5 - f3ii){x cos $ + y sin $ - f3j)) , 
Rl{K^,<!),Y) = ((z-At)(ycos$-xsin$)). (3.26) 
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The explicit $- dependence denoted here is a purely geometrical consequence 
of rotating the from the direction of b to the direction of K±. In ad- 

dition, there is an implicit ^-dependence of the space-time variances, (x^x^) 
= {x^x^){Kj_,^,Y). This ^-dependence characterizes the dynamical corre- 
lations between the size of the effective emission region ("region of homo- 
geneity") and the azimuthal direction in which particles are emitted. Both 
implicit and explicit ^-dependences are mixed in the harmonic coefficients 

R^jJ = ^/_''^4cos(m$)rf$, (3.27) 

= ^f_^Rl^Hm^)d^. (3.28) 

In models in which elliptic deformations dominate and higher than second 
order harmonic coefficients can be neglected, these coefficients satisfy the 
relations [ |175| 



a,{K^, Y) ^ Rt/ ^ li?^/ ^ -RlJ , (3.29) 
a2{K^, Y) ^ R:/ ^ -R:/ ^ -RlJ . (3.30) 

The anisotropy parameter ai vanishes in the absence of position-momentum 
correlations in the source and thus characterizes dynamical anisotropics. On 
the other hand, 0^2 characterizes the elliptical shape of the emission region. 
A violation of the relations ( |3.29| )-( PT^ ) would rule out a large class of model 



scenarios considered to be consistent with the present knowledge about the 
space-time evolution of the collision process. 

Constraints of the type ( p.29[ / |330| ) on the harmonic coefficients lead to a 



minimal azimuthally sensitive parametrization of the two-particle correlator 



T75| : 



C^M^K) ^ 1 + \{K) C,y^{q,K)CM,K,^R) 

xC2iq,K,ijR), (3.31) 

Csymiq, K) = exp[-i?o,o^ ql - Rs/ Qs 

-Ri/qf -2Roi/qoqi], (3.32) 
Ci{q,K,^pji) = exp[-ai(3g^ + g^) cos(<l>-^/'/j) 

+2ai qoQs sin(<l> - ^jr)] , (3.33) 
C2iq,K,ijR) = exp[-a2(go -gDcos2($-V^^) 

+2a2gogsSin2($-^«)]. (3.34) 
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In addition to the azimuthally symmetric part ( |3.32| ) which coincides with the 
Cartesian parametrization, this Gaussian ansatz involves only two additional, 
azimuthally sensitive fit parameters. 







<a^> 







7 8 



Figure 3.2: The HBT anisotropy parameters (oj) as a function of the pa- 
rameter ^ which characterizes the event-by-event reconstruction uncertainty 
in the orientation of the reaction plane. The parameters (ai), (02) are de- 
termined by fitting ( ^.31| ) to an event sample ( p. 36] ) of correlators whose 
reconstructed reaction planes fluctuate around the true impact parameter 
with the probability distribution ( |3.35D . The value ^ = 2 corresponds to a 
reconstruction uncertainty of approximately 30°. 



The major difficulty in determining the harmonic coefficients ( p.27|/^l2^ ) 
from experiment is that they are expected to be statistically meaningful only 
for relatively large event samples while their extraction is conditional upon 
the eventwise reconstruction of the reaction plane. For sufficiently high event 
multiplicities, the probability distributions W{vi,ijjR) of the experimentally 
determined first harmonic coefficients and reaction plane around the most 
likely, "true" values {vi, ipji) is given by a Gaussian of variance 77^ [ p.l6| , |117| , 



■ exp 



vf + vf — 2viVi cos-^ij' 



(3.35) 



An event sample with oriented reaction plane should thus be compared to a 
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weighted average of the parametrization ( |3.31| ), 

Cf^{q,K) = J Vidvid^RW{vi,^R)C^^{q,K). (3.36) 

This effective correlator depends on the event averaged harmonic coefficient 
Vi and the variance rj only via the ratio ^ = vi/rj which is a direct measure 
of the accuracy for the reaction plane orientation [ |116| , |117| , p.65|| . As can 



be seen from Figure ^T^, for an oriented event sample with the typical 30° 
uncertainty in the eventwise determination of the reaction plane angle ipr 
(which translates into ^ ~ 2 [ |165|| ), more than 80 % (50 %) of the anisotropy 



signals ai (02) survive in the actual experimental measurement. Since the 
width 7] is given by event statistics, it is the possible to reconstruct the true 
values ai^2 from the measure values (01,2)- 

Furthermore calculating C^^{q,K) for ^ = 0, i.e., for an azimuthally 
symmetric event sample of finite impact parameter collisions, one can test 
to what extent the azimuthally symmetric HBT radius parameters extracted 
from a fit to such event samples will pick up contributions from non-zero ai 



and a2- This effect is, however, expected to be small [ 165, 175 



3.2 The Yano-Koonin-Podgoretskii parame- 
trization 



The mass-shell constraint K ■ q = 0, explicitly given in (|3.6| ), allows for 
different choices of three independent relative momenta. The Yano-Koonin- 
Podgoretskii (YKP) parametrization, which assumes an azimuthally sym- 
metric collision region, uses the components q± = \J ql + ql, q^ and qi and 
starts from the Gaussian ansatz [ p.82| , |123| , |79| 



(3.37) 



C{q,K) = l + Xexp[-Rl{K)ql-Rl{K){qf-{q'f) 

- {RliK) + E^^iK)) {q-U{K)f . 

Here, U{K) is a (-fC-dependent) 4-velocity with only a longitudinal spatial 
component, 

U{K) = 7(K)(l,0,0,i;(K)) , (3.38) 
1 



7(K) 



(3.39) 
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The combinations of relative momenta {qf — {q'^)^), {q-U{K)Y and qj^ ap- 
pearing in ( p. 37 ) are scalars under longitudinal boosts, and the three YKP fit 



parameters B?\_{K), Rq{K), and -R|(-fC) are therefore longitudinally boost- 
invariant. In contrast to the Cartesian radius parameters, the values ex- 
tracted for these YKP radius parameters do not depend on the longitudinal 
velocity of the measurement frame. This is advantageous in fitting experi- 
mental data. The fourth YKP parameter is the Yano-Koonin (YK) velocity 
v{K) which, as we will see, is closely related to the velocity of the effective 
particle emitter. The corresponding rapidity 

transforms additively under longitudinal boosts. 

Since the ansatz ( p.37|) uses four Gaussian parameters, it is a complete 



parametrization for azimuthally symmetric collisions. These parameters can 
again be expressed in terms of the space-time variances {x^x^) [ |180|| : 



RliK) = Rl{K) = {f){K), (3.41) 
Rl{K) = A-vC, (3.42) 
Rf,{K) = B-vC, (3.43) 




v{K) = I , (3.44) 



where, with the notational shorthand ^ = x + iy, 



2\ 

A = (3.45) 




C ^ (,-Ae))(K). (3.47) 

In these expressions, (y) = (xy) = since we are dealing with azimuthally 
symmetric sources. The kinematical limit K± — > is not free of subtleties, 
as one may guess by finding in the denominator of the above expressions. 
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Indeed, for K± = the mass-shell constraint ( |3.6| ) reads g° = Piqi, and the 
relative momenta q^, qi and q± on which the YKP ansatz ( 3.37|) is based are 



no longer independent. Hence, strictly speaking, the YKP parametrization 
exists only for K± 7^ 0. In practice this does not limit the applicability since 
the Kj_ 0-limit is well-defined for all YKP parameters. 

Mathematically, the Cartesian and YKP parametrizations are equivalent 
and differ only in the choice of the independent relative momentum compo- 
nents. The Cartesian radius parameters can therefore be expressed in terms 
of the YKP ones [ ^ |8| via 



Rl = Rl, (3.48) 
Rl^ = Rl-Rl = /5i7^ {Rl + v'R^^ , (3.49) 
= (1 - A') R\ + 7' (A - [rI + Rfi , (3.50) 

Rli = P±{-PiRl + 7HPi-v){Rl + Rl)) . (3.51) 

This set of equations provides a useful consistency check for correlation data 
analyzed independently with both the Cartesian and the YKP parametriza- 
tions. To invert them, one has to calculate 



1 



A = ^Rls: (3.52) 

B = Rf - ir-Rli + -^-Rdiff) (3-53) 

p± p± 

C = -^Rli + ^Rl,. (3.54) 



and insert them into (|3.41|) - (|3.44|) . These relations imply in particular that 
the YK velocity v{K) can be calculated from measured Cartesian HBT radii. 
In model studies [ |180|| , it was demonstrated that this velocity follows closely 
the velocity of the Longitudinal Saddle Point System (LSPS) which is the 
longitudinally comoving Lorentz frame at the point of highest particle emis- 
sivity for a given pair momentum K. In this sense the YK velocity can 
be interpreted as the effective source velocity. Note that in the Cartesian 
parametrization the kinematical information associated with the YK veloc- 
ity is contained in the cross-term R^i [ |123| , |144|| . 



While the values extracted for Rq{K) and -R|(l^) are independent of the 
longitudinal velocity of the observer system, their space-time interpretation 
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is not. For their analysis the so-called Yano-Koonin frame, which is pair 
momentum dependent and defined by v{K) = 0, offers itself: in this frame, 
the terms ~ vC in ( ^.421 ) / ( ^.431 ) vanish. For a class of Gaussian model emis- 
sion functions including longitudinal and transverse flow it was shown that 
the spatio-temporal interpretation of the flt parameters is then particularly 
simple [H Fl: 



Rh 



K) 
K) 
K) 



(3.55) 
(3.56) 
(3.57) 



In other words, the three YKP radius parameters give directly the transverse, 
longitudinal and temporal size of the effective source in the rest frame of the 
emitter. Especially the last equation ( |3.571) seems to imply that in the YKP 
parametrization the emission duration (t^) can be accessed directly. This, 
however, is model-dependent: in ( 3.56| ) / (|3.57 ) certain terms were omitted on 
the right hand side which can become large in certain model scenarios. For 
example, opaque source models with strongly surface-dominated emission 

m 

1 



give a leading geometric contribution to R^ 



Ra^- 



{{x^) — (y^)) for opaque sources. 



(3.58) 



Large geometric corrections were also observed for transversely expanding 
sources with a box-shaped transverse density proflle [ |163|| . As will be dis- 



cussed in section |5. 4. 1|, flrst checks indicate that opaque source models cannot 
reproduce the experimental data consistently [ p.60| , |177|| , but they play an 



important role in understanding the range of validity of the approximations 

It can happen [ |161| , p.62|| , especially for sources with (x^ — y"^) < 0, that 
the argument of the square root in ( |3.44| ) becomes negative. In this case 
the YKP parameters are not deflned. For such situations a modified YKP 
parametrization was suggested in [ |161| , |162| | which does not have this poten- 
tial problem. The corresponding modifled YK velocity still follows closely the 
fluid velocity at the point of highest emissivity, i.e. also the modifled YKP 
parametrization allows to determine the effective source velocity. However, 



the interpretation of the modifled parameters R'q , -Rji [ |162|| is less straight 



forward than (|3.56|/p37[) ; in particular the parameter R!^^ is in general not 
dominated by the emission duration (t^). 
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So far, the YKP-parametrization has not been extended to colhsions at 
finite impact parameter. 



3.3 Other Gaussian parametrizations 

A plethora of different Gaussian parametrizations can be found in the hter- 
ature. They belong to either of two different classes. The first class contains 
parametrizations which are equivalent to the ones discussed above. A typical 



example for azimuthally symmetric sources is [ 112 1 



C{q,K) = l + X{K)exp[-Rl{K)ql-Rl{K)ql 

-Rl{K)ql-T\K){q'f] . (3.59) 

It provides a perfectly valid azimuthally symmetric ansatz whose four fit pa- 
rameters are, after insertion of the mass-shell constraint ( |3.6| ), seen to be 
in one-to-one correspondence with the Cartesian or YKP ones. One should 
keep in mind, however, that the suggestive notation T'^{K) does not warrant 
a physical interpretation in terms of a temporal extension; also the Rf{K) do 
not only contain spatial information. The interpretation of these parameters 
has to be established again on the basis of space-time variances. Other equiv- 
alent parametrizations can be found in the literature; the relations between 
the various radius parameters are discussed in Refs. [ |162| , |163 . 



The second class contains incomplete parametrizations: either certain 
terms (e.g. the out-longitudinal cross-term in the Cartesian parametrization) 
are neglected, or the ansatz is dimensionally reduced. The prime example is 
the ginv-paranietrization (gf^^ = — {(fY) 

C{q, K) = l + \{K) exp [-i?L(^) ^L] ■ (3-60) 

Here all the different spatial and temporal informations contained in the 
space-time variances {x^x^) are mixed into one fit parameter Rf^^{K), and 
there is no possibility to unfold them again. Furthermore, low-dimensional 
projections of a correlator which is well- described by a complete three-dimen- 
sional Gaussian parametrization in general deviate from a Gaussian shape. 
This is true in particular for projections on q^^^; in fact, it was repeatedly 
observed that C(ginv) is better described by an exponential or an inverse 
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power of ginv For a space-time interpretation of correlation data such in- 
complete parametrizations are not suitable. It is often argued that limited 
statistics forces one in practice to adopt dimensional reductions in the fit 
parameter space. But even then it is preferable to bin the data in three inde- 
pendent g-components first and to project this three-dimensional histogram 
onto different one-dimensional directions for fitting purposes. The parame- 
ters extracted this way can be compared to suitably averaged versions of the 
HBT radius parameters (|3lO|) - (|3T5|) or (|^)-(|^. 

3.4 Estimating the phase-space density 

As shown by Bertsch [ ^ the correlation function can be used to extract the 
average phase-space density at freeze-out. In the present section we describe 
how this works. 

The phase-space density f{x,p,t) of free-streaming particles at time t is 
obtained by summing up the particles emitted by the source function up to 
this time along the corresponding trajectory: 



Here f3 is the velocity of particles with momentum p. For large times t, / 
is normalized to the total event multiplicity, / d^x d^p f{x,p,t)/{27rY = N. 
According to Liouville's theorem, the spatial average of any power of / is 
time-independent after particle production has ceased (t > t/). This is in 
particular true for the average phase-space density 



This quantity can be obtained from the measured one- and two-particle spec- 




(3.61) 




(3.62) 



tra. To this end one calculates, see ( p.24| ). 



vMVi{p2){c{pi,p2)-i) 



^2(Pl,P2) - Vlipi) Vl{p2) 




(3.63) 
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and integrates over q with the mass shell constraint q ■ K = Q. After substi- 
tuting X ^ X + [3x^ (where /3 = K / K^) one obtains 

j d'qS{q-K) [Pi(pi)Pi(p2)(C(pi,P2)-l)] 

X j dx° S{x + f3x\ K) J dy° S{y + I3y^, y\ K) (3.64) 
Jd^xJ:\x,K). (3.65) 



{2nl 
Ek 



In the last step we used the on-shell approximation ^ Ek and introduced 
the time-integrated emission function 

/oo 
dx^ S{x + l3x°,x°,K) . (3.66) 
-oo 



It is easy to show that 
(27r)3" 



j d^x J:''{x, K) = j d^x r{x, K, t>tf). (3.67) 



Combining this with (|3.62|) and (|3.65|) and using the smoothness approxi- 
mation Vi{pi)Vi(j)2) ~ {Vi{K)) on the l.h.s. of ( p.64| ), the phase-space 
density / can be expressed in terms of observable quantities: 

{f){K)^V,{K) Jd\6{q-K){C{q,K)~l). (3.68) 

Using the Cartesian parametrization of the correlator for zero impact pa- 
rameter collisions as given in section |3.1| , the r.h.s. takes the explicit form [ 



17 



Ml cosh Y 



VMK±,Y) = " R,{K)^Rl{K) Rf{K) - {Rl{K)y .(3.70) 



This expression assumes an intercept A = 1 for the correlator. In reality 
a considerable fraction of the observed pions stems from resonance decays 
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after freeze-out. The longlived resonances affect the intercept parameter A 
(see section ^.3.61), and the corresponding decay pions should not be counted 



in the average pion phase-space density near freeze-out. This can be taken 
into account by substituting in (|3.69|) 



ifKK) 
X{K) 



'KK) {f){K) 



(3.71) 
(3.72) 



where the sum in ( p.72|) runs over the resonance fractions fr{K) of longhved 
resonances contributing to the one-particle spectrum at K. 
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Figure 3.3: Rapidity dependence of the average pion phase-space density 



(|3.71|) . The data points are for Au+Au collisions at 10.8 GeV/A measured by 
E877 at the AGS. The statistical errors are smaller than the symbols, the sys- 
tematic normalization uncertainty is of the order of 30 %. The upper (lower) 
lines were obtained assuming a thermal distribution ( ^.73| ) with a temper- 
ature extracted from the high (low) momentum part of the vr" -spectrum. 
(Figure taken from 



112|. 



A first application of this approach was performed by the E877 experiment 
for Au-|-Au collisions at the AGS [ |112| , [T7| . The extracted average pion 
phase-space density in the forward pair rapidity region 2.7 < Y < 3.3 is 
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shown in Figure p.3| as a function of the pair rapidity, averaged over K±. 
As very forward rapidities are approached, it decreases from ^ 0.2 to ~ 0.1. 
This suggests that the phase-space density is largest in the center of the 
colhsion at mid rapidity (here: Y = 1.57) and decreases towards forward 
and backward rapidities. 

More recently a large set of data, including n-p collisions at 250 GeV, 
S-nucleus collisions at 200 A GeV, Pb+Pb collisions at 158 A GeV and 
Au+Au collisions at 10.8 A GeV, was compiled in Ref. [ The resulting 
average phase-space densities are shown in Figure 3.4| . They indicate a strong 
dependence of the spatially averaged phase-space density on the transverse 
momentum K± but very weak dependence on the size of the collision system 
and on the pion rapidity density dN/dj. This latter aspect can be interpreted 
as evidence for a universal pion freeze-out phase-space density in heavy-ion 
collisions [ 

As a simple test of the thermalization assumption, one has compared [ 
112| , |52| the measured phase-space densities in Figs. p.3| and |3.4| to that of a 



thermal Bose-Einstein equilibrium distribution 

f^^(K^, Y) = -— (3.73) 

^ ^' ^ exp[Mxcosh(r-y,)/T] - 1 ^ ^ 

at pair rapidity Y for a source rapidity y^, extracting the temperature from 
the measured single-particle spectra and two-particle correlations. Dynam- 
ical flow effects alter quantitative aspects of f^^{K), but qualitatively the 
phase-space density thus calculated compares surprisingly well with the ex- 



perimentally determined one, see Figs. 3.3 and 3.4. 
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Figure 3.4: Spatially averaged phase-space density in narrow rapidity win- 
dows as a function of the transverse momentum K±. The heavy- ion data 
span more than an order of magnitude in rapidity density dN/dj, but the 
resulting freeze-out densities show much less variation. The i^_L-dependence 
of (/) can be well parametrized by an exponential function. The dashed lines 
show Eq. ( p.73|) for T = 80, 120, and 180 MeV, respectively. (Figure taken 
from [ |52| . ) 
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Chapter 4 

Beyond the Gaussian 
parametrization 



The space-time variances characterizing the Cartesian HBT radius parame- 
ters ( p.l(]| )-( prT5| ) can be written as second derivatives of the correlator ( |3.5| ) 
at q = 0: 



(4.1) 

q=0 



These curvature terms coincide with the experimentally determined half 
widths of C{q,K) only if the correlators is a Gaussian in q. Realistic two- 
particle correlation functions show, however, more or less significant devi- 
ations from a Gaussian shape. The consequences are two-fold: the corre- 
sponding space-time variances do not agree exactly with the fitted radius 
parameters, and the Gaussian radius parameters do not contain all the in- 
formation contained in C{q,K). Nevertheless, qualtitatively all statements 
made above about mixing of spatial and temporal information in the HBT 
radius parameters remain valid. The reason is that the Fourier exponent 
g ■ a; in ( |2.24|) can be written as q ■ {(3t — x). Hence, the q dependence of 
the correlator tests always the same combinations {Pit — Xi) of spatial and 
temporal aspects of the emission function S{x,K), irrespective of the shape 
of the two-particle correlator. 

If the correlator deviates from a Gaussian shape, one can either seek 
a more detailed characterization of C{q,K) supplementing the Gaussian 
radius parameters by a larger set of characteristic parameters, or one may 
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proceed to reconstruct from the measured correlations directly information 
about the emission function without invoking a particular parametrization. 
In the following we discuss both of these strategies. 



4.1 Imaging methods 

Instead of determining information about the emission function S{x,K), an 
alternative analysis strategy [ ^ aims at determining directly from the mea- 
sured true correlator C{q,K) — 1 the relative source function SK^r) intro- 
duced in section ( |2.2.4|) : 



C{q, K)-l = j d\ Kiq, r) SkM . (4.2) 

For the kernel, one usually chooses K{q,r) = \4>q/2{r)\'^ — 1 where 0q/2 
describes the propagation of a pair, which is created with a center of mass 
separation r and detected with relative momentum q. (f)q/2 can include two- 
particle final state interactions, cf. (|2.102|) and (|2.126|) ; for free particle 
propagation, K{q, r) = cos(g ■ r). 



Equation ([4.2| ) can be inverted uniquely, allowing in principle for an un- 
ambiguous reconstruction of the relative source function S'/i:(r). In practice, 
finite measurement statistics on the correlation function and the strongly 
oscillating nature of the kernel K{q, r) render the inversion problem non- 
trivial. Brown and Danielewicz [ |3^ suggest to parametrize S'/<:(r) in terms 
of a finite number of basis functions gj{r), 

N 

SKir)=Y.SAK)g,ir), (4.3) 



and to determine the coefficients Sj (K) from the data. In applications [ ^ 
|37| some insight has been gained on how the reconstruction can be optimized 
by choosing suitable sets of functions gj. 



4.2 moments 

Rather than obtaining the HBT radius parameters from a Gaussian fit to 
the measured correlation function, a quantitative analysis of C{q, K) can be 
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based on expectation values {{g{q))) of the true correlator C{q,K) 
relative momentum space [ |172| , |173|| : 

Jd'qq.qj [C{q,K)-l] 



Rij{K) 



Jd^q [C{q,K)-l] 



2 
ol 



R-li 



R-ls 

?2 



Rli 



Rf 



V Ri 

det R{K)/7t3 J d^q [C{q, K) 



!]■ 



- 1 in 

(4.4) 

(4.5) 
(4.6) 



Similar expressions exist for the YKP parameters [ |172|| . For a Gaussian 
correlator, the intercept parameter A and the HBT radius parameters R^^ 
obtained via these g-moments coincide with the values extracted from a Gaus- 
sian fit. For non-Gaussian correlators the HBT radius parameters and the 
intercept can be defined via ( [4.4| )-( ^l6D . 

Deviations of the correlator from a Gaussian shape are then quantified 
by higher order g-moments, which in turn can be obtained as derivatives of 
the generating function Z{y, K), 

J d'qe^'^ y [Ciq,K)-l] 

{-idy 



Z{y,K) 
{{(ih ■ ■ ■ QiJ 



lnZiy,K] 



y=o 



(4.7) 
(4.8) 



dvh ■ ■ ■ dyi 

It is easy to see that Z{y, K) coincides with the relative source function 

Ziy,K) = SK{y). (4.9) 

So far, only the uni-directional version of these expressions has been used in 
theoretical and experimental investigations [ Restricting the correlator 
along one of the three Cartesian axes, C{K, g^) = C{K, qi, qj^i=0), the cor- 
responding HBT radius parameter and intercept are defined via the relations 
(we use the same notation as for the three-dimensional g-moments) 

1 



RJ(K) 



XAK) 



o,s,l, 



2 m ' 

Jdq,qnC{K,q,)-l] 
Jdq4C{K,q,)-l] ' 

(i?,(K)/v^) / dq,[C{K,q,)-l]. 



(4.10) 

(4.11) 
(4.12) 
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In the case of deviations of C{q,K) from a Gaussian form, the intercepts 
Aj from uni-directional fits can differ in the different directions. A quanti- 
tative measure for the leading deviation of C{q,K) from a Gaussian shape 
is provided by the properly normalized fourth order moments, the kurtoses [ 

ini,ii7i 

^rK) = JM^^l (4 13) 

The uni-directional kurtosis Ai{K) vanishes if the correlator at pair mo- 
mentum K is of Gaussian shape in the gj-direction. The applications of 
g-moments to experimental data is limited by statistics: for higher order q- 
moments the main contribution to the moment integral comes from larger 
relative momenta, while most experimental information is concentrated at 
small q. 



4.3 Three-particle correlations 



The complete set of g-moments (|4.8| ) provides a full characterization of the 
shape of the two-particle correlator. Alternatively, a complete characteriza- 
tion of this shape could be obtained in a Taylor expansion of C{q, K) around 
q = 0. Extending equation ( [4.1|) to arbitrary orders, the corresponding Tay- 
lor coefficients read 



■ (4.14) 

q=0 



dqi^ dqi2 ■■■dq., 

However, due to finite momentum resolution, no detailed information about 
the curvature of C{q,K) at q = is available from the data, and ( [4.14| ) 



cannot be applied in practice. It illustrates, however, clearly that the two- 
particle correlator C{q,K) contains no information about odd space-time 
variances: due to the reflection symmetry C{q,K) ^ C{K, —q), all odd 
space-time variances in ([4.14|) vanish. 



The true two-particle correlator R2 = C — 1 depends only on the modulus 
Pij of the Fourier transformed emission function but not on its phase (pij: 

R2{h3) = C(p„p,)-1 = ^, (4.15) 

pit Pjj 

p^, e'^- = f d^xS (x, \{p, + pj)) e*(P»-f^)- . (4.16) 
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Only the phase (pij contains information about odd space-time 

m 



variances 



(4.17) 



Here Kij, qij denote the average and relative pair momentum for two particles 
with on-shell momenta pi and pj, respectively, and the space-time variances 
( .. )ij are calculated with respect to the emission function S{x, Kij). Three- 



particle correlations give access to this phase as can be seen from [ pO 



'^3lPl,P2,P3j 



^3(Pl,P2,P3) 



^3(Pl,P2,P3) 



/?2(l,2)i?2(2,3)i?2(3,l) 

2cOs(</)i2 + 023 + 03l) , 
C3(Pl,P2,P3)-i?2(l,2) 

-i?2(2,3)-i?2(3,l)-l. 



(4.18) 
(4.19) 



Here -R3 (ra) is the true (normalized) three-particle correlation function. To 
investigate the space-time information contained in the phase combination 
on the r.h.s. of ( [4. IS] ) in more detail one can expand the emission function 
S'(a;, Kij) in ( [4.16|) around the average momentum of the particle triplet: 



K = UP1+P2 + Ps) = UK12 + K23 + K 



31 y 



Kij = K + liqij + Qjk) , 
Using qi2 + ^23 + Qsi = one finds [ |n, |6 



i=^j^k. 



(4.20) 
(4.21) 



'12 



■ Qu Q23 



23 + 031 



OK'' 



^K^' 



[112 <fl2 (I23 + (I23 123 112] 



d^{Xf,)3 d'^{x^)3 d'^{xx)3 



+ 



dK'^dK^ dK^dKi" 



+ 



-2^12 ^23 (912 + to) {x^,XuXx)3 + 0{q) 



^K^'^K^ 

(4.22) 



Here the averages ( . . . )3 have been calculated with the emission function 
S{x, K) taken at the momentum K of the particle triplet. 

The measurable phase $ depends on the odd space-time variances (x^), 
etc., and on derivatives of the point of highest emissivity (a;)3 with respect 
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to K. These reflect the asymmetries of the source around its center. In the 
Gaussian approximation (|3.1| ) they vanish. These considerations show that 
the true three-particle correlator contains additional information which is not 
accessible via two-particle correlations. 

In practice, however, it is very difficult to extract this information. The 
leading contribution to $ is of second order in the relative momenta gjj, and 



in many reasonable models it even vanishes [ Therefore new informa- 
tion typically enters ^^3(^1,^2,^3) at sixth order in q. The measurement of 
the phase $ is thus very sensitive to an accurate removal of all leading q^- 
dependences by a proper determination and normalization of the two-particle 
correlator. These general arguments are supported by model studies which 
found that flow, resonance decay contributions or source asymmetries leave 



generically small effects on the phase [ ^ . 

On the other hand, it was pointed out that the intercept of the normalized 
true three-particle correlator in ( [4.18| ) may provide a good test for the 



chaoticity of the source. Writing the emission function for a partially coherent 
source as = 5'cha+>S'coh, the intercept A3 of r^ is given in terms of the chaotic 



fraction e{p) of the single-particle spectrum [ ^ 



\,{K) ^ r,{K,K,K) = 2^e{K)- , (4.23) 



3 - 2e{K) 
(2 - e{K) 



In contrast to the intercept X{K) of the two-particle correlator, the intercept 



([4.23|) of the normalized three-particle correlator is not affected by decay 
contributions from long-lived resonances which cancel in the ratio ( |4.18| ) [ 
8|. 



Complete small-g expansions of R2 and which generalize the Gaussian 
parametrization ( |3.5| ) to the case of partially coherent sources and to three- 
particle correlations, improving on earlier results [ [3^, 121 ], can be found in 



[ 0. In the framework of a multidimensional simultaneous analysis of two- 
and three-pion correlations they permit to separately determine the sizes of 
the homogeneity regions of the chaotic and coherent source components as 
well as the distance between their centers. 
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Chapter 5 

Results of model studies 



A completely model-independent reconstruction of the emission function 
S{x, K) from measured correlation data is not possible since, due to the 
mass-shell constraint ( 2.27 ), only certain combinations of spatial and tempo- 
ral source characteristics are measurable. Only the time-integrated relative 
distance distribution in the pair rest frame S'i<:(r) can be determined uniquely 
from momentum correlations. This corresponds to a whole class of emission 
functions S{x,K), not a unique one. 

In practice, the mass-shell constraint is not the only problem for the 
reconstruction of the emission function. The statistical uncertainties of ex- 
perimental data turn even the reconstruction of the relative source function 
SK{r) into a complicated task [ 0, where additional model assumptions are 
employed to obtain a convergent numerical procedure. 

Due to these fundamental and pragmatic problems, the analysis of exper- 
imental correlation data starts from a model of the emission function S{x, K) 
from which one- and two-particle spectra are calculated and compared to the 
data. Here one can follow two alternative approaches: either, one simulates 
directly the kinetic evolution of the reaction zone up to the freeze-out stage, 
calculates the one- and two-particle momentum spectra and compares them 
to the data. This is the approach followed in event generator calculations. 
The space-time information is then extracted from within the model simula- 
tion. Alternatively, one applies the tools presented in the preceding chapters 
by using simple parametrizations of the emission function and adjusting the 
model parameters by a comparison to data. This approach makes explicit 
use of the relation between the measured HBT parameters and the space- 
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time variances of the emission function, trying to obtain a direct space-time 
interpretation of the measurements. In this case the dynamical consistency 
of the extracted space-time structure of the source with the preceding kinetic 
evolution of the reaction zone must be established a posteriori. We will here 
concentrate on the second approach; reviewing a large body of model stud- 
ies we illustrate the extent to which momentum correlations have a generic 
space-time interpretation. The insights gained in these models studies are 
summarized at the end of our review into an analysis strategy which allows 
for a simple determination of the relevant space-time aspects of the source 
from experimental data. This strategy is then applied to recent one- and 
two-particle spectra for negatively charged particles measured by the NA49 
collaboration in Pb-|-Pb collisions at the CERN SPS. 



5.1 A class of model emission functions 

We introduce a class of analytical models for the emission function of a 
relativistic nuclear collision, starting from a simple model and then discussing 
several dynamical and geometrical refinements. Models of this class have 



been used extensively in the literature [ 47, 48, 153, |, 42, 44, 171, 172, 175 



5.1.1 The basic model 

Whatever the true particle phase-space distribution of the collision at freeze- 
out is, we expect that its main characteristics can be quantified by its widths 
in the spatial and temporal directions, a collective dynamical component 
(parametrized by a collective flow field) which determines the strength of 
the position-momentum correlations in the source, and a second, random 
dynamical component in momentum space (parametrized by a temperature). 

A parametrization which is sufficiently flexible to incorporate these fea- 
tures but still allows for an intuitive physical interpretation of its model 
parameters assumes local thermalization prior to freeze-out at temperature 
T and incorporates collective expansion in the longitudinal and transverse 
directions via a hydrodynamic flow field u^{x). The source has a finite geo- 
metrical size in the spatial and temporal directions, encoded in transverse and 
longitudinal Gaussian widths R and A77 as well as in a finite particle emis- 
sion duration Ar. Here r = \/t'^ — denotes the longitudinal proper time 
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and T] = [(t + z) / {t — z)] the space-time rapidity. The parametrization is 
optimized for sources with strong, approximately boost-invariant longitudi- 
nal expansion for which freeze-out occurs close to a hypersurface of constant 
longitudinal proper time r = tq. It is thus more suitable for high than for 
low energy collisions. The source is defined by an emission function for each 
particle species r [ |153|, g |2|, ||, |71|, |17| : 



2Jr + l 

(27r)37r Ar 



X exp 



■ m_|_ cosh(y — 77) exp 



p ■ u[x) — fir 



T 



7] 



r 



(5.1) 



2i?2 2{Ar]y 2(Ar)2 

For explicit calculations, it is helpful to express the particle four-momentum 
p^ using the momentum rapidity y and the transverse mass m± = Jp\ + ni^- 
This allows for a simple expression of the Boltzmann factor in ( |5.1| ), 



p ■ u{x) 
T 



(m_L coshy,p_L, 0, m_L sinhy) , (5.2) 

^ cosh(y - r]) coshr^t - sinh?7t . (5.3) 
1 I r 



For sharp freeze-out of the particles from the thermalized fluid along a hy- 
persurface S(x) = {tq cosh 7], x,y,To smb. rj), we would have to choose the 
emission function proportional to p ■ n{x), where 



n^(x) 
p ■ n{x) 



d^a^{x')6^''\x-x') 



m_L cosh(y — t]) 6{r — tq) 



(5.4) 
(5.5) 



The four- vector n^{x) points normal to the freeze-out hypersurface. The 
term m_|_ cosh(y — rj) in the emission function ( p.l|) stems from this geomet- 
rical condition, while we have replaced the 5-function in ( |5.5| ) by a properly 
normalized Gaussian to allow for a finite emission duration At. The factor 
2Jr + 1 accounts for the spin degeneracy of the emitted particle, and a chem- 
ical potential fir allows for separate normalization of all particle yields. The 
ansatz implies that all particles are assumed to freeze out with the same geo- 
metric characteristics and the same collective flow, superimposed by random 
thermal motion with the same temperature. 

For the flow profile we assume Bjorken scaling [ ^ in the longitudinal 
direction, vi = z/t; this identifies the flow rapidity r^gow = |logY 



^ with 
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the space-time rapidity. Assuming a linear transverse flow rapidity profile of 
strength r// in the transverse direction, the normalized flow field u^{x) reads 

u^{x) = (^cosh 7] cosh T]t, ^ sinh r)t, ^ sinh r]t, sinh 77 cosh rjt^ , 

Vt{r)=Vf^- (5.6) 

In spite of the longitudinal boost-invariance of the flow, the source as a whole 
is not boost-invariant unless the longitudinal Gaussian width Arj — > cxo. 

The model emission function (|5.1|) is thus completely specified by six 
common and one species-dependent model parameters: 

T ,r]f,R,Ar],AT,To,Hr. (5.7) 

The number of model parameters can be reduced by assuming chemical equi- 
librium at freeze-out. This provides the following constraint between the 
chemical potentials: 

fir = brfJ-B + SrfiS ■ (5.8) 

Here br and Sr are the baryon number and strangeness of resonance r, and fiB, 
Us are the two independent chemical potentials required for baryon number 
and strangeness conservation in the reaction zone. 



5.1.2 Model extensions 

Comparative model studies investigate to what extent geometrical or dynam- 
ical assumptions put into the emission function S{x, K) leave traces in the 
observed one- and two-particle spectra, and how this allows to distinguish 
between different collision scenarios. To this end one compares, for example, 
the model ( ^.1|) with modified model assumptions about the particle produc- 
tion and emission processes in the collision region. Here we focus on a few 
possible extensions of the basic model { f).l\j which all have a clear physical 
motivation and which have been investigated in the literature. 



Opaque sources 

The particle production and freeze-out mechanisms in heavy ion collision are 
largely unknown. In particular, it is not settled whether hadronic freeze- 
out resembles more the surface evaporation of a hot water droplet or the 
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simultaneous bulk freeze-out leading to the decoupling of photons and the 
transition from an opaque to a transparent universe in the early stages of 
Big Bang cosmology. The parametrization ( |5.1| ) reflects the second type of 
scenario, with a bulk transition from opaqueness to transparency at proper 
time Tq±At. Surface emission can be included into ( ^.1|) by multiplying this 
emission function with an exponential absorption factor [ |8^, |161|| 



Sop^que{x,p) = S{x,p) exp 



-^8/^ (/eff/Ainfp) 



(5.9) 



eflf 



_ y2 POO _ ^12 

hs{r,<f)) = e j e ^ dx' , (5.10) 



where y = r sin (p, x = r cos (p. This extra factor suppresses exponentially 
emission from the interior of the emission region. The particle propagates 
in the x (out) direction. The Gaussians in the expression ( |5.10| ) parametrize 
the matter density seen by the particle according to the geometrical source 
distribution in (15.11). 



Temperature gradients 

The model emission function S{x,p) presented in (|5.1| ) assumes that all vol- 
ume elements freeze out at the same temperature T{x) = T. Since freeze-out 
is controlled by a competition between the local expansion and scattering 
rates, and the latter have a very strong temperature dependence, this is not 



an unreasonable assumption [ 146, 108 1. With this assumption all position- 



momentum correlations in the source in ( [5.1| ) stem from the collective dynam- 
ics characterized by the flow u^{x) in the Boltzmann factor. To contrast this 
scenario with models in which the x-i^-correlations have a different origin, 
model extensions with a particular temperature profile have been studied [ 



160 



' ^ '^ + d' ^ .J ■ (5-11) 



T(x) To V 2i?V V 

Hereby transverse and temporal temperature gradients are introduced via 
two additional fit parameters a and d. The model (|5.11|) implies that the 
production of particles with larger is more strongly concentrated near 
the symmetry axis and average freeze-out time. 
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Emission functions for non-central collisions 



The emission function ( p.l| ) shows ay ^ —y symmetry in the os/-coordinate 
system. It describes an azimuthally symmetric emission region which is ade- 
quate for zero impact parameter collisions. To investigate the new qualitative 
features introduced by collisions at finite impact parameter, one can study 
azimuthally asymmetric extensions [ ITS, ^ of the model emission func- 
tion (|5.1| ). Here we take x and y to be defined in the laboratory system, with 
X in the reaction plane. An elliptic geometric deformation of the source in the 
transverse plane, characterized by an anisotropy parameter e^, is obtained by 
replacing in 



exp 



2R\ 

= Ry/T 



exp 



X 



2PI 
Py 



(5.12) 
(5.13) 



An elliptic deformation of the transverse fiow pattern can be introduced by 



Uf,[X) 



(7_L cosh ?7, Ma-, Uy, 7_L sinh r/) 



7± 



l + ul + ul 



(5.14) 
(5.15) 
(5.16) 



These modifications implement in a simple way some aspects of finite impact 
parameter collisions in the mid-rapidity region. They do not encode for the 
fact that in the fragmentation region the particle emission is peaked away 
from the beam axis. Hence, the total angular momentum L of the system. 



Li = eijk{{xjPk)) = eijk J J d'^xXjPkS{x,p) 



(5.17) 



vanishes for the prescriptions ( |5.12| )-( [5.16| ) given above. A simple parametri- 
zation of the emission function in the fragmentation region which leads at 
least to a finite expression ( p. 171 ) , reads [ |175 



Px 

Py 



x + XY 
R 

XJCOScp) y/T 



(R 

{R + XJ cos (f) y/TTTs 



(5.18) 
(5.19) 
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These dynamical ( p.l8| ) and geometrical ( |5.19| ) assumptions effectively shift 



the center of the particle emission in the transverse plane as a function of 
the longitudinal particle rapidity y with some asymmetry strength x- They 
break explicitly the 180° rotation symmetry of the emission function Slab in 
the transverse plane, which is left unbroken by (|5.12|) - (|5.16|) . 



5.1.3 Resonance decay contributions 

A significant fraction of the most abundant candidates for interferometric 
studies, charged pions, are produced by the decay of unstable resonances 
after freeze-out. To a lesser extent the same problem exists also for kaons. 
Longlived resonances can escape to quite some distance from the original 
freeze-out region before decaying. They then lead to HBT radius parameters 
which are larger than the width of the particle production region. Further- 
more, due to the resonance decay phase-space, secondary pions populate 
mainly the low momentum region and can thus introduce an additional pair 
momentum dependence of the two-pion correlator. To obtain realistic esti- 
mates for the geometry and dynamics of the particle emitting source, reso- 
nance decay contributions therefore must be analyzed quantitatively. 

In the following discussion we focus on charged pions. We include all rel- 
evant resonance decay channels r in the model emission function by writing [ 

61, |33|, Pi 



S^{x,p) = Sf\x,p) + ^ Sr^^{x,p) . (5.20) 

The emission functions 5'r^7r(x, p) for the decay pions are calculated from the 
direct emission functions Sf^^{X, P) for the resonances by taking into account 
the correct decay kinematics for two- and three-body decays. Capital letters 
denote variables associated with the parent resonance, while lowercase letters 
denote pion variables. In particular, and $ here denote the transverse 
mass and azimuthal direction of the parent resonance, in contrast to the rest 
of the review where M± and $ are associated with the pair momentum K±. 
We follow the treatment in [ |157| , |33| , p.72|| . The resonance r is emitted 



with momentum P at space-time point X'^ and decays after a proper time 
r at a;'^ = + if i'^to a pion of momentum p and {n — 1) other decay 
products: 

r — > n + C2 + C3 + ... + Cn- (5-21) 
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The decay rate at proper time r is Fe"^"^ where F is the total decay width of 
r. Assuming unpolarized resonances with isotropic decay in their rest frame, 
Sr^Tr{x,p) is given in terms of the direct emission function S^^^{X, P) for the 
resonance r by 



M 



X 



dsg{s) 



6{p- P - E*M) / drVe 



-Ft 



5f (X, P) 



(5.22) 



Variables with a star denote their values in the resonance rest frame, all 
other variables are given in the fixed measurement frame, s = {J27=2Pi)'^ 
is the squared invariant mass of the {n — 1) unobserved decay products in 



( |5.21| ); it can vary between s_ = {Y17= 



and s_ 



(M 



g{s) IS 



the decay phase-space for the (n 
decays it reads 

9{s) 



1) unobserved particles. For two-particle 



Anp* 



-5 [s — 



whereas the three-particle decay phase-space is given by 

Mh 



Q{M, 171,1712,7713) 



[77l2 + 771^ J 



{7712 - 



2718 

ds' 



Q{M, 777, 1712, ms) 



(5.23) 



(5.24) 



L 



m] 



X ^s^ — s'y^S- — s'^J (m2 — m^y — s' . (5.25) 

Equation (|5.22| ) can be simplified considerably: for p± 0, the energy- 
momentum conserving 5-function constrains the angle $ of the resonance 
momentum = (M_l cosh y, Pj_ cos $, Pj_ sin $, M_l sinh y) to one of two 



orientations if its decay product propagates in the out-direction: 

y ^($-$±) 

± P±P± sin $± ' 
7n±M± cosh(y -y) - E*M 



5{p-P- E*M) 



cos $ 



± 



PlPi 



This allows to do the ^-integration in ( [5.22| ), leading to 



Sr-^-K^X-! p) 



(5.26) 
(5.27) 

(5.28) 
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J2± sums over the two azimuthal directions in ( |5.27| ), and /j^ indicates the 
remaining integrations over the resonance momenta; for details see Ref. [ 



5.2 One-particle spectra 

The one-particle momentum spectrum, determined as the space-time integral 
( |1.3|) of the emission function S{x,p), is sensitive to the momentum distri- 
bution in S{x,p) and thus allows to constrain essential parts of the collision 
dynamics. It contains, however, no information about the space-time struc- 
ture of the source. Statistical errors on one-particle data are significantly 
smaller than those on the two-particle spectra. In practice, exploiting the 
temperature, flow and resonance mass dependence of the one-particle spec- 
trum therefore allows to reduce the model parameter space significantly even 
before comparing model predictions to the measured two-particle correla- 
tions. 



5.2.1 Transverse one-particle spectrum 

We begin by discussing the rapidity-integrated transverse momentum spec- 
trum of the models discussed in the last subsection. Both direct "thermal" 
pions produced in the collision region and those stemming from resonance 
decays contribute to the spectrum, see (|5.28|) : 



dm?, dm?, , 

For the model this expression takes a compact form [ |157] , [145| , |1 72|| : 

dN^'' 2Jr + 1 



(27ri?2 . 2toAt]) e'^'-/^ d ^ 



dMl 47r2 

X (^coshr7i(0) h (^sinhr7i(0) • (5.30) 

^ = r/R is the rescaled transverse radius. Obviously the geometric parame- 
ters R, Arj, To of the source enter only in the normalization of the spectrum. 
The product of the spatial and temporal extensions of the thermal source 
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determines the total particle yield, but not its momentum dependence. To 
further constrain these parameters is not possible without using two-particle 
correlations. According to ( p.29[ / p75(]| ), the m^-dependence of the pion spec- 
trum is fully determined by the temperature T (or T(^) if T is r-dependent), 
the rest masses Mr and chemical potentials fir of the resonances, and the 
transverse flow profile r]t{^) = rjf^"'- 

To study the dependence of one- and two-particle spectra on the com- 
position of the resonance gas, the resonance fractions fr can be computed [ 



172 



Mp) 



Jd^xSr^^{x,p) _ dNl/d^p 



(5.31) 



As we discuss in section |5.3.6| , these fractions play an important role in esti- 
mating the correlation strength \{K). In Fig. |5.1| , we plot the pion transverse 
mass spectrum dN.„/dm\ and the resonance fractions fr(j,P±) of the model 
(|5.1|) for two sets of source parameters. All resonance decay contributions 
are shown separately. The resonances uj, r] and r]' contribute with 3-body 
decays whose decay pions are seen to be particularly concentrated at small 
p± . Comparing the cases of vanishing and non- vanishing transverse flow one 
observes the well-known flattening of the transverse mass spectrum by trans- 

i52|, [TT5|, wm, m wm 



verse radial flow 
an effective "blueshifted" temperature [ |102 



The direct pions reflect essentially 



teff 



1 + (A) 
^ 1 - (A) 



(5.32) 



where the average transverse flow velocity (jSt) is directly related to rjf. This 
clearly does not allow to separate thermal from collective motion. Deviations 
from (|5.32|) are seen for the transverse mass spectra of heavier particles at 
low p± and have been used to determine the temperature and flow velocity 
separately [ 13(:, 158|. One of the main goals of two-particle interferometry 



is to obtain a more direct measure of the transverse expansion velocity at 
freeze-out. 



86 




Figure 5.1: Left: The single-pion transverse mass spectrum for the model 
(|5.1| ) at T = 150 MeV, transverse flow ?7/ = (upper panel) or rjf = 0.3 (lower 
panel), and vanishing chemical potentials Hr = 0. The overall normalization 
is arbitrary, the relative normalizations of the various resonance contributions 
are fixed by the assumption of thermal and chemical equilibrium. Right: The 
resonance fractions fr{j,P±) according to Eq. (|5.31| ) for the same parameters. 
Left column: as a function of transverse momentum at central rapidity. 
Right column: fr as function of rapidity at pj^ = 0. 

5.2.2 Rapidity distribution 

The rapidity distribution is given by 



dy 



J m±dm±d(f) j d'^x S{x;p±, (f),y) . (5.33) 



For the model (|5.1|) its width is dominated by the longitudinal width At] of 
the source. This is a consequence of the assumed boost-invariant longitudinal 
flow profile. Since the resonance decay fractions fr are essentially rapidity 
independent (see Fig. p.l| ), resonance decay contributions do not significantly 
affect the rapidity spectrum. 
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5.2.3 Azimuthal dependence 

For collisions with non-zero impact parameter the triple-differential one- 
particle spectrum ( p.24| ) contains information about the orientation of the re- 
action plane. The harmonic coefficients Vn which characterize this azimuthal 



dependence are given in terms of the Fourier transforms [ 165, 166 



fln, bn) = r2n AN ' i^'^^) 



/o^-E£ (cos(n0), sin(r20)) # 

r2n 7-, dN 

h 

an = Vn cos{nipR) , bn = Vn sin(n?/'/j) . (5.35) 

According to ( p.34| ) they are normalized to the azimuthally averaged double 
differential particle distribution, and Vq = 1. A symmetry argument similar 
to that employed in ( |3.2CI| / P^^TD implies that in the limit p± — > the 0- 



dependent terms vanish, the emission probabilities in different azimuthal 
directions become equal, and 

lim Vn(p±) = for all n > 1. (5.36) 

Px->-0 

Furthermore, the odd harmonic coefficients f „ vanish at midrapidity for sym- 
metric collision systems, due to the remaining — > + tt symmetry in 
the transverse plane. A more explicit calculation depends on details of the 
model. For the emission function ( |5.12| )-( |5.16| ) one can show that V2 oc p\ 



for small values of p±. For small transverse flow, the leading dependence on 
the anisotropy parameters and is given by 



This describes correctly the main features of a numerical study of this model 



175| |. Both geometric and dynamical deformations manifest themselves in 
the single-particle spectrum only for expanding sources with rjf 0. The 
relative minus sign in the numerator of ( |5.37| ) reflects the different signs in the 



definitions of and e/ in ( |5.13|) and ( |5.15| ). Once this is taken into account, 



( |5.37 ) shows that an increasing spatial deformation eg or an increasing flow 
anisotropy ef lead to similar effects on the azimuthal particle distribution. 
Therefore they cannot be separated without also using information from two- 
particle correlations. 



88 



5.3 Two-particle correlator 



In this section we study the question how characteristics of the emission 
function are reflected in the momentum-dependence of the measured particle 
correlations. This was discussed already in section |^ in the context of the 
model-independent relations between the space-time variances of S{x, K) 
and the HBT radius parameters of C{q,K). According to this discussion 
spatial and temporal geometric information about the source is contained in 
the q-dependence of the correlator, while the pair momentum dependence 
characterizes dynamical properties. These statements can be made more 
explicit in the context of specific model studies. 

For quantitatively reliable studies of "realistic" (i.e. sufficiently complex) 
emission functions the determination of the two-particle correlator requires 
a numerical evaluation of the Fourier integral in ( |1.4|) . On the other hand, 
simple analytical approximations for the HBT radius parameters allow to 
summarize the main physical dependencies of the measurable quantities in 
an intuitive form and are quite useful for a qualitative understanding. We 
give a combined discussion of both the analytical approximations and the 
exact numerical results for the model (|5.1| ). 



5.3.1 Saddle point approximation of HBT radius pa- 
rameters 



Characterizing the emission function S{x, K) by the Gaussian widths 
{B~^)^^ (K) = {Xf^Xu){K), is more generally apphcable than a saddle point 
approximation around its center x{K) which was earlier suggested [ |153| , ^ 



The tensor B^y{K) characterizes essential features of the phase-space sup- 
port of the emission function even if S{x, K) is not different iable or if its 
curvature at the saddle point does not represent its average support suffi- 
ciently well. In all these cases the Gaussian widths B^y{K) still translate 
directly into Gaussian radius parameters, as discussed in section |^. Never- 
theless, a saddle point approximation of S{x, K) can be technically useful 
for the evaluation of the integrals ( p.3|) when approximating the HBT radius 
parameters as averages over the emission function [ |153|, ^ 4S , 171 . 



To illustrate the use and limitations of the analytical expressions thus ob- 
tained one can study a simplification of the model (|5.1|) which is particularly 
amenable to analytical calculations: we neglect the resonance contributions 
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and set Arj = oo and At = 0. This represents a source with exact longitu- 
dinal boost-invariance and a sharp freeze-out at time tq [ 153| , |171|| . For 



such a source the cross-term Roi vanishes in the LCMS (see ( |3.17| )); it will 
be discussed in the following subsection. 

Expanding the exponent of the emission function S{x, K) given in (|5.1| ) 
around = one obtains in this case 

In the saddle point approximation the terms bilinear in y specify the (inverse 
of the) "side" radius parameter Rg, while the "outward" radius parameter 
receives an additional contribution from the finite emission duration (t^) = 
ro(sinh^r7) according to [ |^ 

RliK±) = R'M±) + myf3lr'o. (5.39) 
Rli^^) = TT^' (5-40) 

These simple expressions illustrate several of the key concepts employed in 
HBT interferometry: the overall size of the transverse radius parameters is 
determined by the transverse Gaussian widths of the collision region, and the 
difference R^ — R^ is proportional to the emission duration j3j_{P). (Even 
a sharp freeze-out at r = tq corresponds to a finite region in t = tq cosh r) 
since the source distribution is sampled over a finite longitudinal range.) 
Most importantly, however, the radius parameters R^ and R^ are sensitive 
to the transverse flow strength i]f of the source: the HBT radius shrinks 
for finite rjf since a dynamically expanding source viewed through a filter of 
wavelength K is seen only partially. This shrinking effect increases for larger 
values of M± proportionally to the ratio The M^-dependence of Rg 

is a consequence of transverse position-momentum correlations in the source 
which here originate from the transverse collective flow. 

Saddle point integration also leads to simple expressions for the longitu- 
dinal radius parameters. Due to boost-invariance, Ri is /^/-independent, and 
one has to calculate Rf = {z^) = Tq ((sinhr/^) — (sinhr^)^). We summarize 
the results for different approximation schemes: 

^ [to O (-^) , [106]] (5.41) 



R2 ~ 



Ml 
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R 




to O 



rp2 



m 



(5.42) 



[exact for r// = 0, [ H] (5.43) 



In general the longitudinal radius is proportional to the average freeze-out 
time To and falls off strongly with increasing transverse momentum. In 
this case the transverse momentum dependence signals longitudinal position- 
momentum correlations in the source. Compared to the transverse directions, 
it is much stronger since the source expands predominantly in the beam di- 
rection. 

According to the Makhlin-Sinyukov formula ( p. 411 ), fitting an 1/M±- 
hyperbola to (extracting the temperature e.g. from the one-particle spec- 
trum), the source parameter tq can be determined. Early estimates of tq have 
been obtained by this argument. The improved calculations ( p.42|) and ( p.43| ) 
show that for realistic temperatures of the order of the pion mass corrections 
cannot be neglected. Similar remarks apply to the slopes of the "side" and 
"out" radius parameters ( ^.39|/|5^ ). Especially, all the results presented 
here are based on an expansion around = 0, while the true saddle point 
Xfj, for finite ?7/ is shifted in the out-direction; this can be seen from (|5.38|) . 
In Ref. [ |171|| an approximation scheme was developed which takes these 
saddle point shifts into account and allows to systematically derive improved 
expressions for the HBT radii (|5.39|) -( ^.42|) . The resulting expressions are 
involved and a numerical evaluation of the radius parameters is often more 
convenient. 

The sources for quantitative uncertainties of the saddle point approxi- 
mation are two-fold: The determination of the saddle point = {x^) is 
done only approximately in many model studies. The resulting inaccuracies 
turn out to be substantial and put severe limitations on the quantitative 
applicability of the analytical expressions quoted above. This is illustrated 
in Fig. Moreover, whenever the Gaussian widths of the correlator devi- 
ate significantly from its curvature at q = 0, the relation between space-time 
variances and HBT-radii becomes quantitatively unreliable, see chapter ^. In 
the following we therefore use the above simple analytical expressions only for 
qualitative guidance, basing a quantitative discussion on numerical results. 

Equations (|5.39|) - (|5.43|) have been used to extract the transverse flow pa- 
rameter rjf and the "freeze-out time" tq. This is not without danger. As 
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numerical calculation of 
space-time variances 

Gaussian width of nume- 
rically calculated correlator 



result of saddle point 
approximation around"x=o. 



K^(MeV) 



Figure 5.2: i^'x-dependence of the HBT radius parameters Ro (a), Rg (b) and 
Ri (c) for the emission function ( |5.1|) with T = 150 MeV, tq = 3 fm/c, R = 3 
fm, At] = oo, At = 0. Curves for different values of the transverse flow rif 
are shown. Differences between numerically evaluated space-time variances 
and Gaussian widths indicate deviations of C{q, K) from a Gaussian shape. 
Differences to the dash-dotted lines reflect the limited validity of a naive 
saddle point approximation around = 0. 



pointed out in Ref. [ in this particular model the M^-dependence of the 
HBT radius parameters reflects the longitudinal and transverse flow velocity 
gradients. These are responsible for the reduced homogeneity lengths com- 
pared to the total source size. The parameter rj^ in ( |5.4(]| ) really reflects the 
transverse flow velocity gradient [ 42|, and the parameter Tq in ( ^.41|) -( ^.43 ) 



similarly arises from the longitudinal velocity gradient at freeze-out. From 
the general considerations of section |3.1| we know that the absolute position of 
the source in space and time cannot be measured. Therefore, strictly speak- 
ing, To cannot be directly associated with the absolute freeze-out time. Such 
an interpretation of tq makes the additional dynamical assumption that the 
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longitudinally boost-invariant velocity profile existed not only at the point 
of freeze-out but throughout the dynamical evolution of the reaction zone. If 
this were not the case and the collision region underwent, for example, longi- 
tudinal acceleration before freeze-out, the real time interval between impact 
and freeze-out would be longer. 



5.3.2 The out-longitudinal cross-term 

The cross-term R^i of the Cartesian parametrization vanishes in the LCMS 
for longitudinally boost-invariant systems or in symmetric collisions at mid- 
rapidity. This follows from the corresponding space-time variance ( p.l3| ) 
which vanishes in the LCMS if the source is symmetric under z ^ —z. This 
reflection symmetry is broken, however, in the forward and backward rapidity 
regions for the systems with finite longitudinal extension. 




-0.15 -0.1 -0.05 0.05 0.1 0.15 2 * 0.6 0.8 1 



Q,„„JGeV) (GeV/c2) 



Figure 5.3: Left: Contour plot of the correlation function in the qo-qi-p\.a.ne 
obtained by NA35 in S+Au collisions at 200 A GeV. The data are for forward 
rapidity pairs, 3.5 < Y < 4.5, the momentum differences are evaluated in 
a system at rapidity 3. (Figure taken from [|^.) Right: M^-dependence of 
the out-longitudinal radius parameter for the model ( |5.1|) at = 1-5 in 
the LCMS. The curves are shown for different values of the transverse flow, 
rjf = 0.4 (solid) and rjf = 1.0 (dashed). (Figure taken from [ |163|| .) 



As first observed by NA35 [0], following a proposal of Chapman et al. [ 41 



J^, the Gaussian correlator is then an ellipsoid in q-space whose main axes 



do not coincide with the Cartesian ones: it is tilted in the out-longitudinal 
plane (see the l.h.s. of Fig. |5.3|). This tilt is parametrized by the size and 
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sign of R^i- In the forward rapidity region, R^i is positive in the LCMS and 
negative in the CMS [43, 13, 147]. The sign is reversed for the backward 
rapidity region. 

Also, R^i has a characteristic transverse momentum dependence. At van- 
ishing Kj_, the X —X symmetry of the source is restored and R^i vanishes, 
see ( p.21|) and the Lh.s. of Fig. ^.3| . At large K_i_, on the other hand, the 
homogeneity region peaks sharply around the point of highest emissivity, and 
the z —z symmetry is again approximately restored [ |180|| . The resulting 



generic M± dependence is seen in Fig. |5.3| : 
reaches a maximum and then decreases again. 



rises sharply at small Kj 



5.3.3 The Yano-Koonin velocity 

The Yano-Koonin-Podgoretskii parametrization [ |182| , p.23| , |180|| describes 
the correlator of an azimuthally symmetric source by a longitudinal velocity 
v{K) and three Gaussian radius parameters R^lK), R'^^{K), and Rl{K). 
The Yano-Koonin velocity v{K) contains important information about the 
longitudinal expansion of the source. To investigate this, a detailed study [ 
18CI|| was done within the framework of the class of models (|5.1| ), of the 



relation among the following different longitudinal reference frames: 

• CMS: The centre of mass frame of the fireball, specified by r^o = 0. 

• LCMS (Longitudinally CoMoving System [Q): a pair-dependent fra- 
me, specified hj jSi = Y = 0. In this frame, only the transverse velocity 
component of the particle pair is non- vanishing. 

• LSPS (Longitudinal Saddle-Point System [ @): The longitudinally 
moving rest frame of the point of maximal emissivity for a given pair 
momentum. In general, the velocity of this frame depends on the mo- 
mentum of the emitted particle pair. For symmetric sources the point of 
maximal emissivity ( "saddle point" ) coincides with the "source centre" 
x{K) defined in ( |3.2| ). In this approximation, for a source like ( |5.1| ), 
the LSPS velocity is given by the longitudinal component of u'^{x{K)). 

• YK (Yano-Koonin frame [|7^): The frame for which the YKP velocity 
parameter vanishes, v{K) = 0. Again, this frame is in general pair 
momentum dependent. 
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The velocities (or rapidities) of the CMS and LCMS frames can be deter- 
mined experimentally, the first from the peak in the single particle rapidity 
distribution, the second from the longitudinal momentum of the measured 
pion pair. The velocity of the LSPS is determined by the flow velocity at 
x{K) and hence, it is not directly measurable (neither the one- nor the two- 
particle spectra depend on x{K)). However, the YK- and LSPS-systems co- 
incide as long as the particle emission is symmetric around x{K) [ |123 , 18C . 
In model studies based on the parametrization (|5.1|) asymmetries are found 
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Figure 5.4: Calculation of the Yano-Koonin rapidity for the model (|5^ 
without resonances. Source parameters: i? = 3 fm, Tq = 3 fm/c, Ar = 1 
fm/c, A?7 = 1.2 and T = 140 MeV. (a) YK rapidity as a function of the 
pion pair rapidity Y (both measured in the CMS frame of the source), for 
various values of the transverse pair momentum Kx. and for two values of the 
transverse flow rapidity r^j. (b) Same as (a), but shown as a function of 
for different values of Y . The curves for negative Y are obtained by reflection 
along the abscissa, (c) The difference Y^^ — ^lsps between the rapidity of 
the YK frame and the longitudinal rest system of the saddle point, plotted 
in the same way as (a), (d) Same as (c), but shown as a function of for 
different values of Y . 



to be small [ |18CI||. This is seen in the difference between the corresponding 
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rapidities, Y^^^ — I^lsps' plotted in Figure |5^,d as a function of K± and Y 
respectively. The difference is generally small, especially for large transverse 
momenta [ |18(J|| where thermal smearing can be neglected: 



v{K) ^ v,,,,{K) . (5.44) 

The observable YK velocity thus tracks the unobservable LSPS velocity. This 
is important since longitudinal expansion of the source leads to a character- 
istic dependence of the LSPS velocity ^Lgpg on the pair rapidity which - 
using ( ^.44|) - one can confirm by measuring the YK-velocity. Two extreme 
examples illustrate this: for a static source without position-momentum cor- 
relations the rapidity of the LSPS is independent of the pair rapidity Y and 
identical to the rapidity of the CMS: 

■^Lsps ~ o T~~ — ~ const. for a static source. (5.45) 

t — "^LSPS 



In contrast, for a longitudinally boost-invariant model (|5.1| ) with Ar] = oo, 
the longitudinal saddle point lies at t] = Y. In this case, the LSPS and the 
LCMS coincide, 

y^sps = V = ^ for a longitudinal boost-invariant source. (5.46) 

For the model ( |5.1D the measurable YK rapidity plotted in Fig- 

ure |5.4]a,b as a function of the longitudinal pair rapidity Y and transverse 



momentum K±. The plot confirms a linear relation 

Fyj^ const. X Y, for the model (5.1). (5.47) 

with a proportionality constant which approaches unity for large Kj_. Since 
Yyj^ ^ ^LSPS, this linear relation between the rapidity Y^k of the Yano-Koonin 
frame and the pion pair rapidity F is a direct reflection of the longitudinal 
expansion flow. The linear relation shown in Fig. |5.4| a has been confirmed 
subsequently by experiment [ O. There are also first experimental hints for 



the ii'_L-dependence shown in |5.4] a. 

We emphasize that the observation of a linear relation Y^^^ = Y cannot 
be generally interpreted as evidence for boost-invariant longitudinal expan- 



sion, as was suggested in Refs. [ |7|, |180| , In Fig. ^]4|a one sees that 



for the model ( p.l| ) this linear relation is the better satisfied the larger the 
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transverse momentum of the particle pair. This can be easily understood 
as a consequence of the reduced thermal smearing at large Kj^. Whenever 
the random component in the momentum distribution becomes small, the 
observed pair velocities track directly the flow velocity of the emitting vol- 
ume element, irrespective of the actual flow velocity profile. The generic be- 
haviour shown in Fig. |5.4| a therefore indicates longitudinal expansion which 
is sufficiently strong to overcome the thermal smearing, but not necessarily 
boost-invariant. 



5.3.4 Yano-Koonin-Podgoretskii radius parameters 



We discussed in section |33 that, as long as certain asymmetries of the source 
are negligible, the YKP-radius parameters have a particularly simple space- 
time interpretation. The transverse radius parameter always coincides with 
the source width in the side-direction, R\{K) = (y^), and [ |7^, |180 



RlliK) = [[^-f^3:j j - ^{y') ^ {z') , (5.48) 
RliK) = lit- ) - ^if) ^ (?) . (5.49) 



The approximation in these equations results from dropping terms pro- 
portional to (zx), (xt), and (x^— y^). The first two of these terms vanish if the 



source is symmetric around its point of highest emissivity x{K). In Figure 5^ 



the effective emission duration y {t'^) and longitudinal size of homogeneity 
(z^) are compared to Rq and R\\ for the model ( ^.1|) . The approximations 



( |5.48|/p749D become exact for vanishing transverse fiow while differences oc- 
cur especially in Rq for large K± or significant transverse fiow. These can be 
traced to the terms —2{xi) /P± + {x^ — y^) / Pj_ which are neglected in ( |5.49| ). 



For the model ( ^.1|) the space-time variances (xt) and (x^ — y'^) indeed van- 
ish for K± where the azimuthal x-?/-symmetry of the source is restored. 
However, when divided by powers of (3±, the corresponding terms in ( p.49| ) 
result in small but finite contributions even for K± — > 0. For opaque sources 
( |5.9| ) the term —2{xi)/i3± + {x^ — can be the dominant contribution; 

this can lead to large negative values of R^{Kj^). The approximation ( ^.49| ) 
thus breaks down for such opaque sources [ ^ |16(]| , |161|| , and the leading 
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Figure 5.5: YKP radii for the same model as in Fig. |5]J. (a) Rq and 



(i^), evaluated in the YK frame, as a function of M± for three values of 
the transverse flow rapidity rjf, for pion pairs with CMS rapidity y = 0. (b) 
Same as (a), but for pions with CMS rapidity Y = 3. (c) and (d): Same as 
(a) and (b), but for _R|| and the longitudinal length of homogeneity {P), 



evaluated in the YK frame. For Y = 0, and y {z^) agree exactly because 
(3i = in the YK frame. 



i^_L- dependence can be recast in the approximate expression (|3.58|) . The ex- 
perimental data exclude large negative values for Rl{K±) and thus rule out 
certain opaque emission functions [ |160|, |177], |161| . 



5.3.5 Azimuthal dependence of HBT radius parame- 
ters 

Both the calculation of HBT-radius parameters in the saddle-point approxi- 
mation and the numerical calculation have been extended \ 117511 to the finite 



impact parameter model described in (|5.12|) -( p.l6|) . For this model the first 
harmonics vanish, and the zeroth and second harmonics can be written in the 
saddle point approximation in terms of an average size R and a dimensionless 
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anisotropy parameter a2- 
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Figure 5.6: First (upper panel, -R^^i) and second (lower panel, -R^^2) har- 
monic coefficients of the transverse HBT radius parameters of the model 
(|5A8|)/(p39D and the model (p7[4D/(p7[5D, respectively. The subscript "*" 



in -R^i, -R^2 stands for the out (dash-dotted), side (dashed) and out-side 
(dotted) components. Solid thin and thick lines denote the zeroth harmonics 
of the out and side radius parameters, respectively. All calculations are for 
T = 150 MeV, To = 5 fm/c, Ar = 1 fm/c, Ar^ = 1.2 and i? = 5 fm. 



Rsfi = , (5.50) 

o,2 - -^^s,2 - --"-05,2^ = ■ (5.51) 



pc 2 
— JrC„ 



a2 is related to the parameter a2 in ( p.30| ) by a2 = a2R^- According to ( p.51j ) 
the relations ( p.30| ) are exact in the saddle point approximation, with 



R' 



a2 



R^ 1 + ^r^} (1 - el 
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(5.52) 
(5.53) 
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For vanishing anisotropy = e/ = the side radius parameter in ( p.52| ) 
reduces to the corresponding expression ( ^.401) . For non-zero or e/ the 



average size R depends only weakly on these anisotropics. 

In Figure we show the results of a numerical study of the model 
( |5.14|/pTT5D which qualitatively confirms the results of the saddle-point ap- 



proximation. The deviations from the leading dependence 2 '• -^s 2^ ■ 
—R^os 2^ = 1 : — 1 : — 1 of ( p. 511 ) are seen to be small. Similarly, the numeri- 
cal study of the model (|5.18|/prT9D which includes directed transverse flow in 



the forward rapidity region allows to confirm the relation ( |3.29 ) between the 
first harmonic coefficients, -Rq 1^ : R'^s^x '■ Rls,i^ = 3 : 1 : — 1. 

5.3.6 Resonance decay contributions 

Resonance decay pions affect the two-particle correlator by reducing its inter- 
cept X{K) and by changing its q-dependence. Typical examples are shown 
in Fig. p.7| . The modifications of C{q, K) due to resonance decays are a con- 



sequence of the exponential decay law in ( ^.28|) , which provides the emission 
function S'^^tt with a non-Gaussian tail in coordinate space. The latter is 
reflected in a non-Gaussian shape of the correlator. The details depend on 



the lifetime of the corresponding parent resonances [ |172 



Short-lived resonances, F > 30 MeV: In the rest frame of the particle 
emitting fluid element these resonances decay very close to their pro- 
duction point, especially if they are heavy and have only small thermal 
velocities. This means that the emission function Sr^n of the daughter 
pions has a very similar spatial structure as that of the parent reso- 
nance, S^", although at a shifted momentum and shifted in time by 
the lifetime of the resonance. Since the Fourier transform of the direct 
emission function is rather Gaussian and the decay pions from short- 
lived resonances appear close to the emission point of the parent, they 
maintain the Gaussian features of the correlator. 

Long-lived resonances, F ^ 1 MeV: These are mainly the r) and rj', 
with lifetimes cr ~ 17.000 and 1000 fm, respectively, and the weak 
decays of Kg and the hyperons which on average propagate several cm. 
Even with thermal velocities these particles travel far outside the direct 
emission region before decaying, generating a daughter pion emission 
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Figure 5.7: Two-pion correlations for the model (5J.) without transverse 



flow {rjf = 0), calculated according to ( p.20| ). Curves show the correlator 
without resonance contributions (thin solid lines), including pions from p 
decays (long-dashed), other shortlived resonances A, K*, S* (short-dashed), 
from the uj (dash-dotted), and including pions from the longlived resonances 
r], 1]', Kg, S, A (thick sohd lines). 



function Sr^jr with a very large spatial support. The Fourier transform 
Sr~*TT{q, K) thus decays very rapidly for q 0, giving no contribution 
in the experimentally accessible region g > 1 MeV. (This lower limit in 
q arises from the finite two-track resolution in the experiments.) The 
decay pions do, however, contribute to the single particle spectrum 
Sr^iriq = in the denominator and thus "dilute" the correlation. 

In this way long-lived resonances decrease the correlation strength A 
without affecting the shape of the correlator where it can be measured. 

• Moderately long-lived resonances, 1 MeV < F < 30 MeV. There is 
only one such resonance, the u meson. It is not sufficiently long-lived 
to escape detection in the correlator, and thus it does not affect the 
intercept parameter A. Its lifetime is, however, long enough to cause a 
long exponential tail in S^^t^^x, K). This seriously distorts the shape 
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of the correlator and destroys its Gaussian form. 



The main effects of resonance decay contributions on the A intercept param- 



eter are mimicked by the "core-halo-model" [ ^] which assumes that the 



emission function can be written as a sum of two contributions 

S{x,K) = S,{x,K) + Sh{x,K) . (5.54) 

Here, the "halo" source function Sh is regarded as the sum over the longlived 
resonance contributions which are wide enough to be unresolvable by HBT 
measurements. Sh thus affects only the intercept parameter. The "core" 
emission function describes the contributions from the direct pions and short- 
lived resonance decay pions which are emitted from the same central region. 
The model ( |5.54| ) is thus a simplified version of ( |5.20|) and provides a simple 
qualitative picture for the intercept: 

X{K) ^(l- Y: friK)] . (5.55) 

y r= longlived J 

The core-halo model neglects contributions from moderately longlived reso- 
nances, essentially the u, for which the distinction into core and halo does 
not apply. Due to their non-Gaussian shape these affect the fit parameters 
for the intercept A considerably and can lead to quantitative corrections of 
( |5.55|) of up to 10 % [ |172|] . An important consequence following already 



from ( 5.55 ) is the transverse flow dependence of the intercept parameter. 



depicted in Figure p^ . With increasing transverse flow, the i^'x-dependence 



of the A-parameter becomes flatter [ p.64|] , though it does not vanish (see 



Fig. |5.8|). This is important since current measurements are consistent with 



a X^-independent intercept parameter [ |T^, |147 . 

Aside from this (model-independent) lifetime effect, which generically 
increases the effective pion emission region, various model-dependent features 
can affect the degree to which resonance decay contributions change the shape 
of the correlator. In the model ( ^.1| ) for example, the size of the effective 
emission region in the transverse plane shrinks with increasing transverse 
flow rjf approximately according to ( 5.40|) . The simplified core-halo model [ 



49| neglects the proper resonance decay kinematics and thus does not describe 
this effect. As a consequence, the direct resonance emission function 5*^'^ in 
?1]) has a smaller effective emission region than that of the thermal pions. 



102 




200 



400 



800 

Mj^- in MeV 



Figure 5.8: The intercept parameter as a function of M±_ for different 
transverse flow strengths rjf. We use the model ( [5 .11 ) with temperature 
T = Ib^MeV . The calculation is not based on a proper fit of the two-particle 
correlation function, but on ( 5.55 ), including in the sum over resonances the 
contributions from a;, r/, r]\ and K^. 



due to the larger transverse mass M±_ of the parent resonance ("transverse 
flow effect"). For finite transverse flow this reduces the size of the resonance 
emission region and counteracts the lifetime effect. Other models of heavy 
ion collisions [ |141| , |143|| do not show this behaviour and lead to significantly 
different ^-dependences of the correlator. 



5.3.7 Kurtosis of the correlator 

If, as in the presence of resonance decays, the correlator deviates from a 
Gaussian shape, the characterization of C{q, K) via HBT radius parameters 
is not unambiguous. Fit results then depend on the relative momentum 
region covered by the data, on the statistical weights of the different q-bins, 
and on details of the fitting procedure. 

The g-moments discussed in section allow to quantify deviations of the 
correlator from a Gaussian shape [ |173| , |98[| . They are sensitive to differences 
between model scenarios which cannot be distinguished on the basis of Gaus- 



sian radius parameters. To illustrate this point we show in Figure 5.9 both 
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the side radius parameter Rs{K±) and the corresponding one-dimensional 
kurtosis As{K±), calculated according to ( |4.10|) and (|4.13|) . In the present 



case, the side radius parameter calculated from the inverted second g-moment 
coincides rather accurately with the radius parameters extracted from a fit 
to 

C{qi,K) = 1 + Xe-^''^' , i = o,s,l (5.56) 

in the range qi < 100 MeV. For a transverse flow between r]f = and r]f = 0.3, 
Rs shows approximately the same K±-slope. Hence, once resonance decays 
are taken into account, scenarios with and without transverse flow cannot be 



distinguished unambiguously on the basis of (|5.40|) . The physical origin of 
the K±-s\ope of Rg is, however, different in the two situations, and this shows 
up in the kurtosis of the correlator. Without transverse flow, resonance decay 
contributions increase Rg due to the lifetime effect. For non-zero transverse 
flow, on the other hand, the K±-slope arises from the i^'^-dependent shrink- 
ing ( ^.40|) of the effective transverse emission region which is more prominent 
for resonances than for thermal pions. S^" is spatially more extended in the 
transverse plane than S"^"^, and thus "covers" a substantial part of the expo- 
nential tails of Sr^tt [ |1 72|| . As a consequence, the total emission function 



( |5.20 ) shows much smaller deviations from a Gaussian shape for the scenario 
with transverse flow and results in a more Gaussian correlator. This explains 
why the kurtosis As{K±) plotted in Fig. |5.9| provides a clearcut distinction 
between the two scenarios. 



5.4 Analysis strategies for reconstructing the 
source in heavy-ion coUisions 

A realistic emission function of heavy-ion collisions should simultaneously 
reproduce the spectra and particle yields of all observable particle species. 
Every model implies certain constraints between all these observables. For 



instance, in the model (|5.1| ) the assumption of resonance production with 
thermal abundances inside the same space-time geometry relates the pro- 
duction of different particle species. While this may be sufficient to model 
gross properties of particle production, extensions involving additional model 
parameters may be needed to account for finer details (for example, partial 



strangeness saturation or rapidity dependent chemical potentials [ |155|| ). In 
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Figure 5.9: The inverted uni-directional second g- variance Rg of ( [4.10| ) and 
the kurtosis of ( 4.13| ) as a function of K± at mid rapidity for the model 
( |OD with T = 150 MeV, R = 5 im, Ar] = 1.2, tq = 5 fm/c, Ar = 1 fm/c 
and vanishing chemical potentials. Left: rjf = (no transverse flow). Right: 
rjf = 0.3. The difference between the dashed and solid curves is entirely 
dominated by cu-decays. 



this review, we will stay with a simple model and try to describe only the 
rough features of the freeze-out process. We will restrict our discussion to the 
single-particle spectra of negatively charged particles and to two-pion corre- 
lations, setting all chemical potentials to zero. We aim to extract from the 
data the phase-space properties of the pion production region, characterized 
by the model parameters in (|5.1| ). 

5.4.1 Determining the model parameters of analytical 
emission functions 

The model parameters of an analytical emission function should be deter- 
mined by a multi-parameter fit of the corresponding one- and two-particle 
spectra ( |1.3| )-(pr5D to the data. For sufficiently complicated emission func- 
tions, where no accurate analytical approximations of (|1.3|) - (|1.5|) are avail- 
able, this is a significant numerical task. In a model like ( |5.1|) , however, 
certain model parameters are almost exclusively determined by particular 
properties of the measured particle spectra. Based on this observation we 
outline here a strategy for the comparison of ( |5.1|) to data, which uses heavily 
the model studies presented in sections |^]^ and f).'d\ It allows to determine 
the model parameters T, rjf, R, At], At, tq by a significantly simpler method, 
according to the following steps [ p.77|| : 
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Figure 5.10: LHS: x contour plot of a fit to the NA49 /i~-spectrum [ ^ 
Dashed hues are for constant values oirj'j/T. RHS: different combinations of 
temperature T and transverse flow rjf can account for the same one-particle 
slope. 



1. The transverse single-pion spectrum dN/dm\ determines the blue-shif- 
ted effective temperature Teg. 

For pions with T m^, the slope of dN/dm\ is essentially given by 



Teff = Ty'(l + - (A)), as argued in section |^:^ and Fig. ^ 

Resonance decay contributions affect the local slope of dN/dm\ and 



thereby the fit parameter [ 145 |. They have to be properly taken into 
account. Then a fit to the corresponding one-particle spectrum deter- 
mines a "valley" of parameter pairs T and rjf (related to (A)) all of 
which can account for the same data. This is clearly seen in the x^-plot 
of a fit to recent NA49 /i~-spectra [ 0, presented in Fig. p.lQ . 



2. Combining the single-particle spectrum dN/dm\ and the transverse 
HBT radius parameter Rs{Mj_) = R±{Mj_) disentangles temperature 
T and transverse flow rjf. then fixes the transverse extension R. 
In the saddle-point approximation ( p.4U| ) the M^-slope of the transverse 
radius R± is proportional to ^ij/T. In Fig. |5.10| we have superimposed 



lines of constant Tfj/T onto the fit results for dN/dm\. Due to the 
different correlation between T and rjf, the additional information pro- 
vided by -R_l(Mx) allows to disentangle temperature T and transverse 
flow effects [ p.47|| . Once the transverse flow is fixed, the overall size of 
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R±{M±) determines the Gaussian width R of the source according to 
(|5.40|) . Quantitative details change if the saddle-point approximation 
is abandoned, but the qualitative argument survives. In a numerical 



calculation including resonance decay contributions [ |177|| one extracts 
then from the combination of R±{M_i_) and dN/dm\ the following val- 
ues for the model parameters, see Fig. |5.12| : rjf ^ 0.35, T ^ 130 MeV, 
i? 7 fm. 

The single-particle rapidity distribution dN/dj fixes the longitudinal 
source extension Arj. 

The single-particle rapidity distribution ( |5.33| ) is determined by Arj, 
the only parameter which breaks the longitudinal boost-invariance of 
the source (pTl|). Fig. |5.11| shows that the h~ rapidity spectrum is a 
Gaussian with a width of 1.4 rapidity units. In the parametrization of 
the present model, this translates into Ar] = 1.2. 

Pb + Pb, NA49 Preliminary 



• Measured 


h" (n", K", p) 


O Reflected 






CT = 1.4 




<h'> = 680 







Figure 5.11: The rapidity distribution of primary negative hadrons. The 



preliminary NA49 Pb+Pb data are taken from [ ^ . 



4. Ru determines tq. 



In principle, and Ri depends on tq, Arj and Ar [ |171|| . Model 
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calculations [ |161|| show, however, that for ( |5.1| ) the dependence on 



the emission duration Ar is weak. As argued above, can be fixed 
from the single-particle rapidity distribution. The data presented in 
Fig. ^.12| b then clearly favour a value of tq ~ 9 fm/c. In this plot Ar 
was chosen to Ar = 1.5 fm/c. From the arguments given at the end 



of section |5.3.1| , this value of tq is likely to provide a lower estimate for 
the total lifetime of the collision region. 

5. Rq discards opaque sources. 



For the model ( |5.1| ) the YKP-parameter Rq is mainly sensitive to the 
mean emission duration of the source, since for this model the approx- 
imation ( p. 571 ) is satisfied. The large statistical uncertainties of the 
NA49 data for Rq do not allow to constrain the model parameter space 
further, see Fig. ^.12| c. Certain models of opaque sources, however, 
which include the opacity factor ( |5.1CI| ), lead according to ( |3.58| ) to a 
negative radius parameter Rq, and can be excluded 
by the present data. 
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Figure 5.12: Yano-Koonin-Podgoretskii HBT-radius parameters. The fig- 
ure uses preliminary NA49 Pb+Pb data for h~^h^ (squares) and h^h'^ (dia- 
monds) correlations [ |T^ . Final data have since been published in [ |1^ . 

5.4.2 Uncertainties in the reconstruction program 

We now list some sources of uncertainties in the reconstruction program 
described above: 
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1. Particle Identification: 

The above analysis is based on /i^ (all negative liadrons) and (all 
positive liadrons) spectra and correlations. These are dominated by the 
corresponding charged pion contributions. The effect of other particle 
species (e.g. kaons or (anti)protons) on the one-particle spectra can 
be included in ( |5.29| ) and has a negligible effect on the shape of the fit 



(which is what matters) shown in Fig. |5]^. For the two-particle cor- 
relations misidentified particles lead to the counting of "wrong pairs" 
which do not show Bose-Einstein correlations. This decreases the in- 
tercept parameter whose value was, however, not used in the recon- 
struction program. There are no indications that the problem of parti- 
cle misidentification affects the i^^-dependence of the HBT parameters 
significantly. Nevertheless, it would be preferable to use spectra of iden- 
tified particles. This would additionally allow to compare different ef- 
fective source sizes e.g. from pion-pion and kaon-kaon correlations and 
thereby investigate the question whether all particle species are emitted 
from the same source volue. Such identified two-particle correlations 



were measured at the AGS [ |112|| and by the NA44 Collaboration at 
the CERN SPS [|6^. Both of these data sets are, however, restricted to 
narrow regions in K± and Y which limits their usefulness for the above 
reconstruction program. If the different acceptance of the two experi- 
ments is properly taken into account, the vr" 7r~-correlations from NA44 
and the h~ /^"-correlations from NA49 are compatible within statistical 
errors (P. Seyboth and J. P. Sullivan, private communication). 

2. Coulomb Corrections: 

Coulomb corrections affect the size of the HBT radius parameters as 



well as their i^^-dependence [ |13|, |147| , P8[ . For example, the difference 
i?o — in NA35 correlation data was found to be non-zero only after 



changing from a naive Gamow correction based on (|2.103|) to an effec- 



tive correction ( ^.104] ) which takes the finite source size into account 



0. Furthermore, in the NA49 experiment a proper treatment of the 
Coulomb correction proved essential for a successful check of the con- 
sistency relations ( p.48D -( p.54D between the Cartesian and YKP HBT 
parameters [ |13| . 
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3. Uncertainties affecting the Kj^- dependence of HBT radius parameters: 
Particle misidentification and Coulomb correction are not the only pos- 
sible sources of uncertainty which may affect the momentum slope of 
HBT radius parameters. Other sources on the theoretical side are for 



example multiparticle symmetrization effects [ p.76|| , the dependence 
on a central Coulomb charge [ [1^, the mixing of anisotropy effects 
into the HBT radius parameters resulting from collisions with non- 



zero impact parameter in the central collision sample [ |165| , |175 



or 

multiparticle final state interactions. On the experimental side, resid- 
ual correlations in the mixed event sample used to normalize the two- 
particle correlator or the shape of the experimental acceptance in the 
/r_L-y-plane can affect the presented i^^-dependence of HBT-radius 
parameters. 

One usually argues that these uncertainties are small for the size and K^- 
dependence of the HBT-radius parameters. However, already small differ- 
ences in the transverse slope of R±, for example, affect significantly the opti- 
mal combination of fit parameters (T, rjf) in Fig. |5.1CI| a. We therefore consider 



the systematical error in the slope of R±{K±) to be the most important un- 



certainty in the reconstruction program of section |5.4.1 



5.4.3 Dynamical interpretation of model parameters 

Pions, as most other hadrons, rescatter during the expansion stage of a heavy- 
ion collision. Their phase-space distribution S'^(x,p) characterizes the geo- 
metrical and dynamical properties of the final freeze-out stage after their last 
strong interaction. The emission function does not contain direct information 
about the hot and dense earlier stages of the collision. However, the emission 
function, as reconstructed from the spectra and correlation data, provides an 
experimentally justified starting point for a dynamical extrapolation back 
towards the earlier stages. To illustrate this point, we discuss here the values 
of the geometrical and dynamical parameters, extracted for the model (|5.1| ) 
from preliminary NA49 data: 

R ^ 7fm, (5.57) 
T ^ 130 MeV, (5.58) 
rif ^ 0.35, (5.59) 
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To ^ 9fm/c, (5.60) 
Ari ^ 1.3, (5.61) 
At ^ 1.5fm/c. (5.62) 

To obtain from these data a dynamical picture of the coUision process, we 
compare first the two-dimensional rms width obtained from the transverse 
width 7 fm, 



r™ = ^(£2 + f) = V2R^10fm, (5.63) 

with the two-dimensional rms widths of a cold lead nucleus. The nuclear 
hard sphere radius i?hs = r^A^^^ with tq = 1.2 fm for lead is R^^ = 7.1 fm. 
The corresponding two-dimensional transverse rms width 



^cold Pb 



^{x^ + = yS/S i?hs ~ 4.5 fm . (5.64) 



From this we conclude that during the collision the system has expanded by 
a factor ^ 2 from the transverse size of the overlapping cold lead nuclei to 
the transverse extension at freeze-out. For not quite central collisions the 
initial nuclear overlap region is in fact expected to be somewhat smaller than 
given by ( |5.64| ). With an average transverse flow velocity of about 0.35 c 
matter can travel over 4 fm in a time of 9 fm/c. This is barely enough to 
explain the observed expansion. This indicates that indeed, as argued before, 
the parameter tq is a lower estimate of the total duration of the collision. 
The rather low thermal freeze-out temperature of 130 MeV (other analyses 
indicate even lower values [ p.63|| ) differs significantly from the chemical 
freeze-out temperature (~ 170 MeV) needed to describe the observed particle 
ratios in these thermal models [ consistent with a long expansion stage. 
A first attempt to extrapolate the final state characterized by ( ^.57|) -( pl62D 



all the way to the beginning of the transverse expansion [ |ST| has led to an 
estimate for the average energy density at this point of e ~ 2.5 GeV/fm^. 



Ill 



Chapter 6 
Summary 



In this work we reviewed the underlying concepts, calculational techniques 
and phenomenological uses of Hanbury Brown/Twiss particle interferometry 
for relativistic heavy-ion collisions. Compared to the astrophysical appli- 
cations of HBT interferometry, its use for relativistic nuclear collisions is 
substantially complicated by (i) the time-dependence and short lifetime of 
the particle emitting source, (ii) the position-momentum gradients in the 
source resulting from the strong dynamical expansion of the collision region, 
and (iii) other dynamical origins of particle momentum correlations which 
have to be subtracted properly to make a space-time interpretation of the 
measured correlation data possible. 

The Wigner phase-space density ("emission function") S{x,K), inter- 
preted as the probability that a particle with momentum K is emitted from 
a space-time point x in the collision region, provides the appropriate start- 
ing point for the analysis of measured HBT correlations from relativistic 
heavy ion collisions. It accounts for the time-dependence and the position- 
momentum gradients of the source. Moreover, other contributions to the 
momentum correlations between pairs of identical particles, for example from 
final state Coulomb interactions, multiparticle symmetrization effects, or res- 
onance decay contributions, can be calculated once the emission function is 
given. In chapter |^ we discussed this in detail, after deriving the basic relation 
(|1.4|) between the phase-space emission function S{x, K) and the measured 
two-particle momentum correlation C(q, K). 

The key to a geometric and dynamical understanding of the measured 
two-particle correlations are the model-independent relations between the 
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space-time variances (Gaussian widths) of the emission function and the HBT 
radius parameters which are extracted from Gaussian fits (|1.5| ) to the two- 
particle correlator. These were derived in chapter |^. In general, the HBT 
radius parameters I^j{K) do not measure the total geometric source size, 
but the size of regions of homogeneity in the source from which most par- 
ticles with momentum K are emitted. These homogeneity regions typically 
decrease in the presence of temperature or flow velocity gradients which lead 
to characteristic position-momentum correlations in the source. Generically, 
the homogeneity regions are smaller for pair with larger transverse momen- 
tum. The -ff-dependence of the two-particle correlator thus gives access to 
dynamical characteristics of the collision region. 

Expressing the HBT radius parameters in terms of space-time variances 
also shows explicitly how spatial and temporal information about the source 
is mixed in the measured momentum correlations. We explained this first for 
the Cartesian parametrization of C{q,K) and its extension to collisions at 
non-zero impact parameter. We then introduced the YKP parametrization, 
an alternative Gaussian parametrization of C{q, K) which is particularly well 
adapted for collision systems with strong longitudinal expansion: (i) Three of 
the four YKP fit parameters are invariant under longitudinal Lorentz boosts, 
i.e. their value do not depend on the observer frame, (ii) The fourth fit pa- 
rameter is the Yano-Koonin velocity which measures the longitudinal velocity 
of the particle emitting source element. Its rapidity dependence allows to de- 
termine the strength of the longitudinal expansion in the collision region, (iii) 
In the particular observer frame in which the Yano-Koonin velocity vanishes, 
the YKP radius parameters for a large class of emission functions cleanly 
separate the longitudinal, transverse and temporal aspects of the source. 
This considerably simplifies the space-time interpretation of the Gaussian fit 
parameters considerably. 

An important aspect of HBT interferometry is that from a combined anal- 
ysis of the single-particle spectra and HBT correlation radii an estimate of 
the average phase-space density of the source at freeze-out can be obtained. 
Applied to heavy-ion data at the AGS and SPS, this method provided evi- 
dence for a universal freeze-out phase-space density for pions. Its transverse 
momentum dependence is in rough agreement with expectations based on 
models assuming thermalization prior to freeze-out. 

The connection between HBT radii and space-time variances of the emis- 
sion function is based on Gaussian parametrizations of the source and the 
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measured particle correlations. Deviations of the two-particle correlator from 
a Gaussian shape can contain additional space-time information which is not 
contained in the HBT radii. In chapter ^ we described refined techniques 
which give access to such additional information. The imaging method and 



the method of g-moments discussed in sections |4.1| and |4.2| characterize in 
different ways the relative source function 5'/<:(r) which is defined as a time 
averaged distribution of relative distances in the source. Identical three par- 
ticle correlations, discussed in section ^73| , give access to odd orders of the 
space-time variances of S{x, K) which drop out in identical two-particle cor- 
relations. 

The main goal of particle interferometry for relativistic heavy-ion colli- 
sions is to extract from the measured momentum spectra as much informa- 
tion as possible about the emission function S{x,K). We explained why a 
completely model-independent reconstruction of S{x, K) is not possible. In 
practice, one therefore must take recourse to a model dependent approach. 
This was illustrated in chapter ^ for a class of analytical emission functions. 
Comprehensive model studies allowed to separate the generic from the more 
model-dependent features. Simple approximate expressions for the HBT ra- 
dius parameters were given which provide a qualitative understanding of the 
dominant geometrical and dynamical effects. Their quantitative accuracy 
was tested by numerical means. 

The size and momentum dependence of the different HBT radius param- 
eters and the transverse momentum slopes of the one-particle spectra were 
shown to depend in general only on one or two of the model parameters of 
the emission function. This allows for a simple analysis strategy for the re- 
construction of the emission function. It was illustrated in section in an 
application to preliminary data from the 158 A GeV lead beam experiment 
NA49 at the CERN SPS. The extracted source parameters are consistent 
with the creation of a highly dynamical system which after impact expanded 
in the transverse direction over a time of at least 9 fm/c with approximately 
one third of the velocity of light before emitting particles at a temperature 
of around 130 MeV. This information about the hadronic emission region 
provides a starting point for a dynamical back extrapolation into the hot 
and dense early stage of the collision, and it can be directly compared with 
the output of numerical event simulations of relativistic heavy-ion collisions. 
Compared to an analysis based on momentum space information only, the 
additional space-time information obtained from particle interferometry pro- 
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vides severe constraints for our understanding of heavy-ion collision dynam- 
ics. 

Remaining uncertainties in the approach were discussed. Better experi- 
mental statistics and further progress in our quantitative understanding of 
effects which influence the pair momentum dependence of the two-particle 
correlator arc expected to lead in the near future to considerable improve- 
ments in the analyzing power of this method. 
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